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PREFACE 


This final report, as issued by the Applications Division of the 
Environmental Research Institute of Michigan (ERIM) under National Aero- 
nautics and Space Administration (NASA) contract NAS3-22892 for the 
Lewis Research Center (LeRC), covers the contract period from September 
1, 1981 through April 30, 1983. The technical representative for the 
contract officer was Mr. Thom Coney of LeRC. The principal investigator 
was Mr. Fred Tanis with important contributions made by Ms. Jacquelyn S. 
Ott, Mr. Thomas Wessling, and Dr. Russell Moll (University of Michigan). 
This research was conducted by the Resources and Technology Department 
of the Applications Division under the direction of Mr. Fred J. Thomson. 

This contract involves developing algorithms to map chlorophyll -a 
and suspended sediment concentrations in Great Lakes waters using the 
Coastal Zone Color Scanner (CZCS). The technical approach is based upon 
the inherent optical characteristics of Great Lakes waters. A number of 
universities and institutions participated in the project. This report 
covers ERIM's activities in the project during the second phase of the 
program. 
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SUMMARY 

A series of experiments were conducted, in the Great Lakes region, 
to evaluate the application of the Nimbus-7 Coastal Zone Color Scanner 
(CZCS) to assessment of trophic state, verification and spatial refine- 
ment of whole lake models, and observation of temporal and spatial dyna- 
mics of phytoplankton. The focus of the present project is the develop- 
ment and testing of atmospheric and water quality algorithms appropriate 
to the Great Lakes as well as evaluation of existing algorithms devel- 
oped for the marine environment. The quantification of chlorophyll and 
suspended sediment in Great Lakes waters by satellite visible radiometry 
is dependent on a thorough understanding of the radiative transfer pro- 
cesses in the atmosphere, at the water surface, and in the water column 
itself. It has been well established that the particulate and dissolved 
substance content of water affects the apparent water color. By sensing 
color with a high signal to noise ratio in narrow spectral bands, CZCS 
provides a means of mapping the water content which has been heretofore 
unavailable from satellite data. Since the air and water radiometric 
effects are coupled, removal of atmospheric effects becomes critical to 
the success of the Great Lakes application. Once the atmospheric path 
radiance has been removed, the radiance which is scattered upward from 
beneath the surface can be observed and used to determine water para- 
meters. The portion of water radiance reaching the satellite from the 
lake surface can amount to as little as two or three percent of the 
total. Further, the variation due to constituent concentration can be 
much less than that due to atmospheric changes. In addition to these 
atmospheric effects, the water problem requires understanding of inher- 
ent optical properties related to measured quantities of chlorophyll -a 
pigment concentration, suspended mineral particles, and dissolved 
organics. 


1 


Previous remote sensing studies of Great Lakes waters employed less 
quantitative approaches and relied heavily on the available surface 
truth measurements [1,2]. Constituent algorithms developed in these 
latter studies were empirical and had little validity beyond the area of 
surface truth, and could not be applied when surface truth were unavail- 
able. In the present study algorithms were sought which can be based on 
optical properties specific to the Great Lakes and which have reduced 
requirements for surface truth. Further they must be able to predict 
surface concentrations of chlorophyll -a and suspended sediments over a 
wide range of values and they must be capable of making predictions over 
water masses which exhibit spatial variation of both atmospheric haze 
and surface concentrations. 

The project was carried out in two phases. The first phase in- 
volved data collection and development of CZCS geometric correction 
software [3]. During the summer of 1980 a series of CZCS overflights of 
the Great Lakes were made with the simultaneous collection of surface 
truth data. The surface truth data consisted of point sampling for water 
quality and optical parameters and airborne flightlines with the Ocean 
Color Scanner. The second phase consisted of data reduction and algor- 
ithm development for both atmospheric corrections and quantification of 
water parameters. ERIM's primary task under this second phase was CZCS 
algorithm development for the Great Lakes. 

The investigation of possible atmospheric correction algorithms 
included evaluation of three existing algorithms and the development of 
a new algorithm based upon the available 1980 CZCS experiments. Because 
of the spatial variations in the concentration of aerosol particles each 
algorithm considered applied a correction on a pixel -by-pixel basis. 
Existing algorithms considered included (1) the Gordon algorithm [4], 
(2) the Smith and Wilson algorithm [5], and (3) the pseudo-optical depth 
algorithm due to Jain [6]. A fourth algorithm developed under this 
project was designated the residual component algorithm. Each of these 


algorithms was adapted to operate in two modes. First, the algorithm 
could be operated in a calibration mode to extract the necessary aerosol 
parameters from a hazy portion of the scene. Second, the algorithm 
could be used with the data derived parameters to process the scene on a 
pixel -by-pixel basis to remove the spatially varying component due to 
aerosol particles. Extraction of aerosol parameters for a given scene 
is dependent on two basic requirements. First, there must be some water 
masses in the scene which have uniformly low and predictable levels of 
chlorophyll-a pigment concentration. It must be possible to reliably 
estimate the upwelling radiance for the calibration areas under a 
variety of illumination conditions. This need can best be satisfied 
with a combination of upwelling radiance measurements or reflectance 
model calculations. Second, the scene must contain sufficient haze over 
one or more of these calibration areas to provide a basis for aerosol 
spectral characterization. 

A comparison study made with these four separate atmospheric algor- 
ithms indicated that the algorithm derived in this project was superior 
to the existing algorithms. Three criteria were used to evaluate the 
performance of each atmospheric correction algorithm. First, derived or. 
measured local optical depths were used with an atmospheric radiance 
model to estimate the aerosol portion of the path radiance. These path 
radiances were compared with those predicted by each algorithm. Second, 
the upwelling CZCS spectra, or that implied from known concentrations of 
principal constituents, were compared with that obtained from each of 
the correction algorithms. Third, the spatial patterns of haze known to 
be present in the raw uncorrected imagery were examined to see if they 
were removed by the algorithm. The residual patterns in the corrected 
imagery were also compared with known or suspected patterns of 
phytoplankton productivity and suspended solids. 

Results from this comparative analysis suggested that the Gordon 
algorithm and its assumption (that the upwelling radiance at 670 nm is 
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zero) are not valid for the Great Lakes waters. Acceptable performance 
could not be obtained even in areas with known minimum productivity. 
Furthermore, it was demonstrated that for extensive areas of the Great 
Lakes the upwelling spectra have a significant component at 670 nm which 
can be utilized to extract information on constituent concentrations. 
Under conditions of light haze each of the three other algorithms 
appeared to be capable of haze removal in CZCS bands one, two, and 
three. The project developed algorithm (4) appears to be capable of 
normalizing the most severe haze conditions. Conditions can exist, 
however, where the haze is sufficiently dense to make recovery of 
surface reflectance impractical. Radiance model studies showed that 
predicted values of the upwelling radiance produced from algorithm (4) 
to be the most valid, but additional surface upwelling radiance data 
would be helpful to further verify water radiance model calculations. 
For the CZCS subscenes tested, algorithm (4) developed under this 
project showed the best overall performance at removing the effects of 
aerosol haze. 

The CZCS water parameter algorithm developed under this project is 
designed to predict both the concentration of suspended mineral parti- 
cles (mg/£) and the concentration of chlorophyl 1-a pigments (ug/t). 
Initially, band ratio and other empirical algorithms developed for 
marine environments were applied to the Great Lakes waters with little 
success [7]. These latter water algorithms were not expected to work 
well for Great Lakes waters since the inherent optical properties of 
these waters are different from those found in an open ocean 
envi ronment. 

Great Lakes water algorithms were developed using a comprehensive 
volume reflectance/atmospheric transfer model. This model is capable of 
predicting CZCS measured radiances at selected viewing geometries with a 
flat water surface and with specified atmospheric and water absorption 
and scattering parameters. This model incorporates Monte Carlo 
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calculations of subsurface reflectance which are related to the 
absorption and scattering properties via a power series expansion [8]. 
The water reflectance model was combined with an atmospheric model so 
that the entire water/atmosphere radiative transfer process could be 
simulated. This model was used to simulate actual CZCS viewing 
geometries and the observed satellite measured radiances in each of the 
CZCS bands. Aerosol and Rayleigh optical thicknesses were estimated for 
each scene date and viewing geometry. Surface truth measurement data 
were used in the water model to estimate the subsurface reflectance. 
The combination of these parameters allowed simulation of the entire 
illuminating and viewing geometry and estimation of the expected CZCS 
radiance. Once the atmospheric parameters were sufficiently 
established, underwater absorption and scattering coefficients could be 
adjusted to improve the fit between model predicted and measured values 
of the CZCS radiance in each band. 

Several different sets of water optical properties were tried based 
on those reported by Jain, Morel and Priur, and Gordon [6,8,9]. 
Underwater optical properties measured in vitro by the NASA Langley 
Research Center during the summer experiments were also incorporated 
[10]. Results from the simulation analysis indicated, as might be 
expected, that no single set of optical properties worked well for all 
waters of the Great Lakes. Optical properties of nearshore waters which 
are often laden with resuspended, eroded, and tributary sediments were 
found to be complex and different between Western Lake Erie, Southern 
Lake Michigan, and Western Lake Superior near Duluth, Minnesota. In 
addition, these waters have different optical characteristics from those 
associated with open waters of central Lakes Michigan, Huron, and 
Superior. 

Surface measurements were to be obtained on all of the Great Lakes 
from other sources, but unfortunately they were insufficient to allow 
derivation of separate water algorithms for coastal and open clear 
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waters of each of the Great Lakes. No surface truth was collected from 
Lakes Huron and Ontario. Under such circumstances a compromise optical 
model for the entire Great Lakes was appropriate and consistent with the 
limited set of coincident surface truth measurements. A direct solution 
search technique described by Hooke and Jeeves [11] was used with the 
simulation model to derive a representative set of optical parameters 
for the Great Lakes. A water parameter extraction algorithm for CZCS 
was based upon this latter optical model. The algorithm consists of a 
series of polynomial prediction equations in the band 1/band 3, and the 
band 1/band 4 ratios which were derived from model simulated radiances. 
Comparisons made with the measurement data showed the algorithm capable 
of predicting concentrations with an average error of approximately 0.5 
Log(C) for chlorophyll-a pigment concentration and for suspended 
sediment concentrations. 

The derived atmospheric and water algorithms were used to extract 
chlorophyll-a concentrations from the available CZCS data for several 
open lake and coastal areas of the Great Lakes. These values were found 
to compare well with those reported in an earlier Upper Great Lakes 
Reference Study which attempted to characterize these waters from 
existing surface measurements [12]. The satellite and measured values 
compared especially well for the open clear waters for which average 
concentrations do not change dramatically from year to year. In produc- 
tive areas such as Saginaw Bay and Western Lake Erie, the CZCS extracted 
values covered a broad range of values which encompassed the measured 
values. For such complex waters it is difficult to characterize the 
trophic state with a single average value without dealing with the 
complex spatial sampling problem. In this regard the CZCS extracted 
values were considered to be potentially more useful for that 
characterization. 

While it is recognized that additional verification and field 
experiments are needed, a great deal of progress has been made toward 



developing quantitative CZCS algorithms for the Great Lakes region and 
similar coastal type waters. Furthermore, it now appears feasible to 
develop stand alone algorithms which once calibrated would require only 
a minimum of surface truth data to have application to multiple scenes 
and location within the Great Lakes. 
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2.0 


INTRODUCTION 

The needs for monitoring and surveillance of the Great Lakes water 
quality have been recognized repeatedly by both Canada and the U.S. 
Federal agencies. This need is mandated by the Great Lakes water 
quality agreement of 1972. Contained within this agreement is the 
requirement for monitoring programs to provide for the assessment of 
water quality conditions and for the effectiveness of corrective 
actions. To meet this requirement, large scale surveillance plans, 
utilizing ship surveys, have been adopted and implemented throughout the 
Great Lakes. Present ship survey techniques are both costly and are 
limited in temporal and spatial resolution. Satellite remote sensing 
offers the potential for increasing both the efficiency and the accuracy 
of shipboard measurements. 

The NIMBUS-7 Coastal Zone Color Scanner (CZCS) was specifically 
designed to quantitatively measure biomass related chi orophyl 1 -pigment 
concentrations in the ocean and, therefore, represents a unique oppor- 
tunity for testing the feasibility and applicability of satellite water 
quality monitoring in the Great Lakes. Previous research and demonstra- 
tion programs have been conducted with aircraft and Landsat have shown 
the ability to measure certain synoptic water quality data and opportun- 
ities to extend ship sampling capabilities [2]. The present approach is 
a quantitative one, whereas previous studies tended to be qualitative, 
involving image interpretation. 

Attempts to extract quantitative water quality data from Landsat 
were based upon empirical models which relied heavily on surface ship 
sampling. The validity of these relationships could not be extended to 
other lakes or satellite scene dates or other water masses within the 
same lake. These difficulties were caused at least by three separate 
sources of error. 
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First, no accounting was made of spatial or temporal variations in 
atmospheric optical parameters. Secondly, the inherent optical proper- 
ties of Great Lakes waters were not measured, modelled or otherwise 
considered in the derivation of remote sensing relationships designed to 
predict constituent concentrations. Thirdly, the spectral bands avail- 
able with Landsat MSS were not ideally placed to make such measurements. 
The present study has addressed each one of these issues in an effort to 
develop extraction algorithms representative to the Great Lakes region. 
Because of the complexity and diversity of Great Lakes water quality 
parameters a credible evaluation for CZCS requires that data be 
collected for the entire range of water conditions for which the CZCS 
data could be applied. The focus of the present program is the 
development and testing of atmospheric and water algorithms appropriate 
to the Great Lakes as well as evaluation of existing algorithms 
developed for the marine environment. 

2.1 STATEMENT OF THE PROBLEM 

The quantification of substances in Great Lakes waters by satellite 
visible radiometry is dependent on a thorough understanding of the radi- 
ative transfer processes in the atmosphere, at the water's surface, and 
in the water column itself. It has been well established that the 
content of water, be it particulate or dissolved substances, affects the 
apparent color. By sensing color with a high signal -to-noise ratio in 
narrow spectral bands CZCS provides a means of looking at the water 
content which has been heretofore unavailable from satellite data. 
Since the air and water effects are coupled to the CZCS radiometric 
data, removal of atmospheric effects becomes critical to the success of 
Great Lakes verification. Once these are removed, the radiance which is 
scattered upward from beneath the surface can be observed clearly by the 
satellite. Effectively, the radiance reaching the satellite from the 
lake surface amounts to a few percent of the total radiation and 
consequently much of the radiation received is from atmospheric 
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backscatter and surface reflectance (see Table 2.1). Furthermore the 
variation in radiance at the satellite due to change in constituent con- 
centration could be on the order of 1% while the corresponding varia- 
tions due to atmospheric changes can be considerably higher. The 
spatially varying atmospheric component is due principally to aerosol 
scattering. The significance of the atmospheric problem for a water 
target has been demonstrated by Hovis and Lung [14] and more recently by 
Quenzel and Kaestner [15] who compared the variability of the atmosphere 
with the reflected light from phytoplankton suspensions. Thus, an essen- 
tial project goal and significant barrier to resolution of quantitative 
information from the water are the reduction of interf erring atmospheric 
effects. 


TABLE 2.1 ATMOSPHERIC INFLUENCE ON 


SIGNALS MEASURED 

OVER THE OCEAN [13] 


Wavelength 

(nanometers) 

Percent Undesired Radiance 
Clear Water Turbid Water 

440 

85.6 

81.9 

520 

82.5 

67.7 

550 

85.5 

65.1 

670 

97.8 

83.6 


Understanding the components of the upwelling radiance which 
contain information on the subsurface constituents requires knowledge of 
both the apparent and inherent optical properties of the water and how 
these properties relate to measured quantities of chi orophyl 1 -a 
pigments, suspended sediments, and dissolved organics. The work of 
Bukata et al [16] on Lake Ontario described the optical properties of 
these components which provided a basic model for the present study. 
Many previous studies relied heavily on the availability of ship 
collected surface truth and employed empirical techniques to relate 
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remote sensing derived variables to surface truth data [17]. Resulting 
algorithms could not be easily applied to other lakes or scene dates. 
For the present study algorithms were sought which could be based on 
optical properties specific to the Great Lakes and which would reduce 
the need for collection of coincident and extensive surface truth. 

2.2 PROJECT PARTICIPANTS 

Prior to beginning, the first phase of the project an ad hoc group 
was organized and coordinated by the NASA Lewis Research Center to 
promote cooperation among the various investigations. The group was 
directed by the NOAA Great Lakes Environmental Research Laboratory 
(GLERL). The principal participants in the project included NASA Lewis 
Research Center (LeRC), ERIM, University of Wisconsin, University of 
Minnesota (Duluth), University of Michigan Center for Great Lakes 
Research, NASA Langley Research Center (LaRC), Naval Oceans Systems 
Center (NOSC), and Science Applications Incorporated (SAI). These 
participants brought a wide variety of backgrounds and capabilities to 
the project. During the second phase of the project ERIM provided 
coordination services to the group by organizing a series of meetings 
and making personal contacts. One of the principal functions of this 
effort was to disseminate the data to all project participants. 

2.3 OBJECTIVES AND APPROACH 

The objectives of the CZCS Great Lakes Experiment were as follows: 

(1) Assure the acquisition of sufficient standardized ground truth 
sampling and ancillary data to fully evaluate the potential of 
CZCS to determine and monitor Great Lakes water quality 
parameters. 

(2) Apply the acquired data pairs to validate CZCS water quality 
algorithms as well as determine their limitations and 
accuracies when used to provide whole lake parameter maps and 
interlake comparisons. 


12 




(3) Apply CZCS products to demonstrate the utility of the data in 
the solution of particular Great Lakes water quality monitoring 
problems. 


The approach taken to meet these objectives consisted of a multi- 
element program consisting of (1) surface sampling, (2) aircraft under- 
flights, (3) algorithm developing including geometric and atmospheric 
correction, and (4) algorithm validation. During 1980 a series of mea- 
surements were made using point sampling from ships and area sampling by 
making a series of aircraft underflights with the ocean color scanner. 
These measurements are discussed in detail in the next section and were 
designated to support the refinement of subsurface volume reflectance 
models. Specific surface sample analysis of interest included chloro- 
phyll and suspended mineral concentrations, spectral signatures, species 
identification, particle size distribution, and transmittance. 

Because large portions of Great Lakes waters exhibit high concen- 
trations of suspended particles, the volume reflectance was expected to 
be significant in the CZCS band 4 (670 nm). It was considered unlikely 
that the NOAA algorithm for atmospheric and surface reflectance correc- 
tions would be accurate for the Great Lakes area. For this reason, 
ancillary data collected on the surface (i.e., atmospheric optical 
thickness) and by aircraft were designed to support and refine an atmos- 
pheric algorithm for the Great Lakes. Correction accuracies were to be 
tested with a number of surface upwelling radiance measurements. 

An essential element to this validation experiment was the precise 
geometric correction of the standard CZCS image data which would allow 
image-to-image registration, and registration of surface truth/satellite 
data pairs. Further, the geometric correction was felt necessary to de- 
velop image products accurately scaled and sized for Great Lakes condi- 
tions and applications. Development of an ephemeris driven scanner 
model was considered key to development of a geometric correction 
algorithm. 
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Atmospheric and bio-optical algorithms developed under this 
program, as well as those developed by the CZCS NET, were to be evalu- 
ated and tested with the available surface truth data pairs. To the 
extent possible it was planned to test the spatial and temporal limita- 
tions of specific algorithms. The applicability of CZCS data for per- 
forming Great Lakes surveillance tasks was to be evaluated by construc- 
tion of suitable product map depicting chi orophyl 1 -a pigment and 
suspended mineral concentrations. These maps were distributed to 
projected participants for evaluation and comment. 

2.4 ENVIRONMENTAL CONDITIONS IN GREAT LAKES WATERS 

Waters of the Great Lakes have distinctive properties from those of 
the deep ocean and to some extent from coastal marine environments. In 
general the deep ocean waters display very low levels of productivity 
with surface chlorophyll-a concentrations which are less than 0.5 yg/fc. 
The optical properties of these waters are dominated by the presence of 
phytoplankton. For these waters, suspended mineral particles are 
essentially not present. These waters have been classed by Morel as 
Type I [9]. Pigment algorithms have been developed by Clark and Gordon 
[7] for these waters which have applicability to large portions of the 
oceans with little need for surface truth calibration. 

The balance of oceanic waters, consisting primarily of coastal 
waters are generally more productive and display values of chlorophyll-a 
pigment which are greater than 0.5 yg/£. These waters exhibit varying 
concentrations of suspended sediments which have a significant and 
sometimes dominant impact on the optical properties. These complicated 
coastal waters have been designated as Type II waters. 

Great Lakes waters exhibit a wide range of optical characteristics 
which are highly variable both temporally and spatially. Surface waters 
of central Lake Superior, northern Lake Huron, and Georgian Bay have 
very low concentrations of suspended mineral particles and have the 
lowest values of chlorophyll-a color pigments. These waters are probably 


of the Type II class but they are very clear and certainly approach in 
optical characterization those classed as Type I. Waters of Green Bay 
of Lake Michigan, Saginaw Bay of Lake Huron, Lake Erie, and much of Lake 
Ontario are highly productive and have as a result high concentrations 
of chlorophyll-a and suspended mineral particles. These waters are 
frequently very turbid and highly reflective, have a shallow optical 
depth, and display surface concentrations which are highly variable both 
spatially and temporally. The balance of Great Lakes waters lie between 
these two extremes. 

Optical characterization is complicated further by the fact that 
the spectral characteristics of suspended mineral particles are differ- 
ent for different areas of the Great Lakes. For example, red silts and 
clays associated with shoreline soils of western Lake Superior near 
Duluth influence the optical properties of off-shore waters and make 
them optically different than those brown and grey sediments found in 
western Lake Erie. Since some of these sediments reach clearer open 
waters of each lake, potentially these waters could also display dif- 
ferent optical properties due to different sediment types. Possibly a 
greater influence in these latter open lake waters is the species 
distribution of phytoplankton. Many of the clear open lake waters tend 
to be dominated by diatoms whereas the nearshore waters are dominated by 
greens and the highly productive areas by blue-green algae. 

Ideally one would want to derive or measure directly the optical 
properties of each sediment type and algal species and perhaps other 
components and create an optically layered model dependent on the 
concentration and distribution of each component. However, the implied 
sea truth requirements for such an analysis are far beyond the scope of 
the present program. Further, the differences in optical properties of 
algal species and suspended sediments have not been thoroughly 
investigated nor has their influence on apparent properties derived from 
remote sensing data been demonstrated. Furthermore, most studies to 



date conducted with Type II waters have shown predicted concentrations 
which are within only a factor of two of measured values. Such results 
are due in part to the decreased sensitivity of the sensor to changes in 
concentration in more turbid and highly reflective waters and to the 
obvious surface sampling problem in highly productive and more turbid 
areas which often display large variations within the spatial resolution 
of CZCS. 

In this context, subtle differences in optical properties will have 
less importance. With a modest amount of surface truth available, as. 
described in the next section, the approach taken in the present program 
was to derive algorithms which would have applicability for the greatest 
extent of Great Lakes waters. 

Because of frequent cloud cover in the Great Lakes area it usually 
is not possible to observe all or even most of the lakes simultaneously 
cloud free, further complicating the collection of adequate surface 
truth for validation, but underlining the need for algorithms which once 
calibrated have a minimum dependence on surface truth. 

In addition to highly complex water optical parameters, the Great 
Lakes region has atmospheric aerosol components which can be expected to 
vary vertically and geographically in both type and concentration. 
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3.0 


DATA COLLECTION METHODS 

A multisensor approach was taken in the collection of surface truth 
data coincident with the collection of CZCS scenes for the Great Lakes. 
This approach involved surface sampling and aircraft underflights with 
the Ocean Color Scanner (OCS). This approach was designed to produce an 
extensive set of optical and biological measurements at a few positions 
in each lake. The aircraft underflights were conducted using an instru- 
ment similar to the CZCS allowing the collection of multiple CZCS-like 
signatures from multiple altitudes over the surface collection sites and 
cross-lake transects which were nearly coincident with the CZCS 
overpass. These latter data were intended to aid interpretation of the 
shipboard point samples and satellite data through variable spatial 
resolution measurements of the spectral distribution of upwelling light 
at both the water surface and at various altitudes. 

3.1 SURFACE SAMPLING 

A brief summary of the Great Lakes data sets collected during 1980 
is presented in Table 3.1. Three principal experiment locations were 
selected: (1) Catawba area of western Lake Erie, (2) Grand Haven area 
in Lake Michigan, and (3) Duluth area of western Lake Superior. Surface 
truth data were collected by the research vessel Hydra (Ohio State 
University) in Lake Erie; the Shenehon (NOAA/GLERL) in Lake Michigan and 
the Leptadora (Univerity of Minnesota) in Lake Superior. Aircraft data 
were acquired by the NASA/LeRC F-106 aircraft instrumented with the OCS. 
Aircraft data were collected at altitudes of 500 and 41,000 feet. The 
500 foot flights were flown over the surface truth vessel in order to 
obtain optical measurements with minimal atmospheric interference. 
Aircraft flights were made within one hour of local noon and usually 
close to the time of CZCS passes over the Great Lakes. Sun glint was a 
special problem for these aircraft data collected near local noon. Sun 
glint interference with the spacecraft sensor was largely avoided by 
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setting the instrument tilting mechanism. Additional aircraft flights 
were conducted for some dates during the morning and afternoon to 
correspond to solar zenith angles between 40 and 60 degrees in order to 
minimize the sun glint problem. The approximate areas of the Great 
Lakes imaged by the aircraft OCS scanner are shown in Figures 3.1 and 
3.2. Figure 3.1 depicts those areas imaged while underflying the Nimbus 
spacecraft. Figure 3.2 shows those areas imaged during the non-satel- 
lite additional flight runs. Aircraft line orientation is dependent on 
the solar zenith angle. The analysis of these aircraft data was part of 
the University of Minnesota study effort and they were not utilized in 
the present phase II study. 




TABLE 3.1 1980 GREAT LAKES EXPERIMENT SUMMARY 

DATA 


No. 

Date 

Surface Sampling 

Flight Altitudes (feet) 

Time 

CZCS Orbit 

1 

11 June 

Lake Michigan 

500, 500, 
500, 41K, 

500, 41K 
41K 

Noon 

PM 

8243 

2 

14 July 

Lake Erie 

41K, 41K, 
41 K, 41K, 

500 

500 

AM 

Noon 

8699 

3 

17 July 

Lake Erie 

500, 41K, 
41K, 41K, 

41K, 41K 
500, 500, 500 

AM 

Noon 

8740 

4 

23 July 

Lake Michigan 

41K, 41K, 

500, 500, 500 

Noon 

8823 

5 

6 August 

Lake Superior 

41K, 41K, 
500, 41K, 

500 

41K, 41K 

AM 

Noon 

9017 


During the time period from July 15 to July 25, 1980 additional 
optical measurements were made at each of the surface collection areas 
by the Naval Oceans Systems Center (NOSC) and NASA/LaRC. The 
measurements included use of the NOSC in-situ submersible radiometer to 
measure upwelling and downwelling irradiance. Measurements were made 
with this radiometer during satellite and aircraft overflights. During 
this same time period the optical parameters, including absorption. 


18 



Superior 





\ 


FIGURE 3.1 AIRCRAF 




FIGURE 3.2. OTHER 



2erjm 

attenuation, and scattering coefficients were measured on twenty-one 
water samples from the three primary experiment sites. Samples were 
also obtained from western, central and eastern Lake Erie as well as 
Green Bay. 

The collection and analysis of all of the bio-chemical data were 
supervised by the University of Michigan Great Lakes Research Center 
Laboratories. Surface truth measurements are summarized in Table 3.2. 

3.2 CZCS DATA COLLECTION 

The CZCS data collected under the 1980 validation experiment 
included five scene dates where substantiating surface truth data were 
collected. These include June 11, July 14, July 17, July 23, and August 
6. The validation analysis conducted for this study included each of 
these five scene dates. Weather conditions for these five dates are 
reported in Table 3.3. Additional CZCS scenes were collected on April 
18, May 22, June 30, July 18, July 25, August 1, September 8, September 
11, and September 18. The June 11 and July 23 scene dates provided the 
best overall cloud free coverage of the Great Lakes. Cloud coverage 
limited the analysis to Lake Erie for July 14 and 17 and to Lake 
Superior for August 6. These latter scenes were not requested by LeRC 
for analysis under the present project. 



TABLE 3.3 WEATHER 

CONDITIONS AT 

PRIMARY SAMPLING 

SITES 

Date 

Visibility/ 
Cloud Cover 

Wind Speed 

Wind Direction 

Wave Height 

June 11 

Clear with haze 

0-3 knots 

180-200° 

< 0.5' 

July 14 

Clear with haze 

7 knots 

225° 

< 0.5' 

July 17 

20% Cumulus clouds 

9 knots 

240° 

1-2’ 

July 23 

Clear with light 
haze 

5 knots 

315° 

1-2* 

August 6 

Clear with light 
haze 

calm 

calm 

< 1.0' 
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TABLE 3.2 SURFACE TRUTH MEASUREMENTS 


Measurements 

Instrument Method 

Sampling Method 

Spectral 

Ranqe 

Downwelling 
Diffuse and 
Global Irradiance 

Tracor Northern 
Spectroradi ometer 

Multiple Deck 
Samples 

400-800 

nm 

Subsurface 
Up/Down Radiance 

NOSC Subsurface 
Multiband Radiometer 

Subsurface 
Measurements at 
0, 0.5, 1.0 secchi 
depths 

440-670 

nm 

Secchi Depth 

30 cm Secchi Disk 




Volume Scattering 
Function 

Modified Brice 
Phenoix Scatterometer 

Large volume 
water sample to 
Secchi depth 

450-800 

nm 

Beam Attenuation 

NASA/LaRC 

Transmissometer 

Same as above 

450-800 

nm 

Absorption 

Coefficient 

NASA/LaRC 
Absorption Meter 

Same as above 

450-800 

nm 

Chlorophyll -a 
and Photo- 
Pigments 

Spectrophotometric 
with acetone 
extraction 

Samples at 0.01, 
0.5, and 1.0 
Secchi depths 
Triple A1 ignots 



Total Suspended 
Sol ids 

Gravimetric 
Ash-free dry weight 

Same as above 



Particul ate 
Organic Carbon 

HP-CHN Analyzer 
Gas Chromatography 

Same as above 



Temperature 

BT-Type/Electronic 

Profile 
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3.3 COLLECTION AND ANALYSIS OF THE NASA LaRC OPTICAL DATA 

The inherent optical properties measured by the NASA/Langley group 
of scientists included absorption, beam attenuation, and the scattering 
phase function. The scattering function was subsequently integrated to 
estimate the total and back scattering coefficients. A few additional 
measurements of this type were made by NASA/LeRC using a separate set of 
instruments which supported those made by the Langley group. The 
measurement instrumentation could not be placed on board the sampling 
vessel so large volume water samples (50 gal.) were transported to the 
NASA/LeRC facilities. For each of the twenty-one optical data sets the 
measured properties pertain to a particular mix of constituents in the 
lake sample. These latter optical data are included in Appendix A. If 
sufficient number of measurements are made in this manner and if the 
principal constituents are known then multivariate regression analysis 
can be used to derive the inherent optical cross sections for a common 
constituent. 

Of the twenty-one optical data sets taken, three Lake Erie samples 
contained sufficient quantities of sediment to saturate the optical 
measurement instrumentation. Samples collected from Green Bay in Lake 
Michigan and from western Lake Superior were found to have distinctive 
local optical properties. Four of the six samples collected from the 
Grand Haven area were essentially sediment free and the presence of low 
concentrations of phytoplankton made absorption measurements difficult 
and the results uncertain. Attempts to include these latter samples in 
the statistical analysis proved to be unproductive. 

The most representative samples, in our opinion, included nine 
samples from western Lake Erie and two samples from Lake Michigan. 
Measured optical properties for these samples are given in Table 3.4. 
Several regression models were formulated and tested against these 
selected optical measurement sets. Regression models were based upon 
chlorophyl 1-a, phaeophytins, total residue, and ashed residue. These 
models are discussed further in section 6.1. 
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^Samples reported here were used to determine optical cross sections 
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3.4 SUBMERSIBLE RADIOMETER DATA 

The NRSC submersible radiometer was used to make upwelling down- 
welling irradiance measurements just above the surface, at the surface, 
and at one and two meters below the surface. Due to wave action the 
subsurface data were found to be absent or inconsistent. The above 
surface radiance measurements and the estimate of the subsurface 
reflectance are compiled in Table 3.5 for Lake Erie sampling stations on 
July 14 and July 17 as well as the Lake Michigan Grand Haven stations on 
July 23, 1980 [18]. Results are reported at 488 and 500 nm in addition 
to the first four CZCS wavelengths. The upwel-ling measurements made on 
July 14 and 17 are generally consistent both in spectral shape and in 
magnitude with model predicted radiances based upon the measured water 
constituents. The Lake Michigan July 23 measurements, on the other 
hand, were found to lack expected spectral shape and were somewhat 
higher than expected. While these data were found useful in the 
verification of optical models and extraction algorithms for the Great 
Lakes there were insufficient measurements for purposes of model 
development. Further comments on these data can be found in the 

discussion of results in section 5.9. 

3.5 DOWNWELLING IRRADIANCE DECK MEASUREMENTS 

A TRACOR Northern spectrophotometer was instrumented on the deck of 
the sampling vessel for purposes of collecting frequent measurements of 
the diffuse and direct solar irradiance. Presurvey calibrations were 
made with an integrating sphere. Preliminary analysis of these data 
indicated that while the measured spectral irradiance was of the 
approximate correct value the shape of the irradiance spectra were peak 
at a longer wavelength than would be expected from model calculations. 
Since measurements were not made precisely at the time of CZCS 
overflight no direct application to the present studies was made. 
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TABLE 3.5 NOSC SUBMERSIBLE RADIOMETER DATA 
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3.6 ATMOSPHERIC OPTICAL THICKNESS MEASUREMENTS 


A NASA/LeRC developed solar radiometer was used to make atmospheric 
total optical thickness measurements. These measurements were made at 
each of the primary sampling sites. The estimate of optical thickness 
was made from a series of solar radiance measurements made over a five 
or six hour period. By measuring only the radiance of the solar disk 
and the solar elevation one can estimate the total optical thickness 
from the following transmittance equation 

E = rE Q e " T/v o (1) 

where t is the total optical thickness, y Q is the cosine of the solar 
zenith angle, E Q is the extraterrestrial irradiance, and r is the 
responsivity of the radiometer. The optical thickness can simply be 
estimated by plotting log E versus l/y Q where the slope becomes the 
estimate of t. Measurements were made at 380, 450, 500, 610, 749, 873, 
and 1040 nm using selected interference filters. Linear interpolation 
was used to estimate the atmospheric optical thickness at CZCS wave- 
lengths. Radiometer optical thicknesses derived by this method are 
shown in Table 3.6 [19]. Subsequent radiometric analysis showed these 
measurements to have the proper spectral shape but incorrect magnitude. 
This result was not surprising since the instruments were not properly 
recalibrated prior to the experiment. 

TABLE 3.6 MEASURED TOTAL OPTICAL THICKNESS 


Date 

Location 

443 nm 

520 nm 

550 nm 

670 nm 

June 11 

Lake Michigan, Grand Haven 

0.605 

0.417 

0.380 

0.284 

July 14 

Lake Erie, Catawba 

0.733 

0.506 

0.449 

0.317 

July 17 

Lake Erie, Catawba 

0.632 

0.412 

0.375 

0.285 

July 23 

Lake Michigan, Grand Haven 

0.441 

0.440 

0.400 

0.284 

August 6 

Lake Superior, Duluth 

0.529 

0.403 

0.373 

0.283 


Note: Values reported here have been estimated by linear interpolation from 

those reported by NASA LeRC [19]. 
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3.7 COLLECTION AND ANALYSIS OF WATER SAMPLES 


Water samples were collected from a series of offshore stations and 
as close as possible to the time of CZCS overpass. The original plan 
was to collect samples from stations located at the apex of a triangular 
array. In actuality the stations were laid out along a line perpendi- 
cular to the shoreline. Figure 3.3 shows the approximate location of 
sampling stations at each of the primary sampling sites. At each 
station samples were collected in triplicate from each of three depths. 
These depths were one meter, one half the Secchi depth, and the Secchi 
depth. When the Secchi depth was less than two meters samples were also 
collected at one third and two thirds this depth. Water samples were 
split into three aliquots and filtered for determination of chlorophyll, 
particulate organic carbon, and suspended solids. Chlorophyll was 
analyzed spectrophotometrically using an acetone extraction procedure. 
Pigments were analyzed for A, A2, B, and C fractions. Suspended solids 
were estimated gravimetrical ly by weighing prepared oven dried filters 
and then ashing the filters to obtain the ash-free dry weight. Parti- 
culate organic carbon samples were analyzed with an HP-CHN analyzer. 
Dry filters were analyzed by infrared gas chromatography for C, and gas 
chromatography for N and H [20]. 

In addition to the primary sampling sites, samples were also 
gathered from the non CZCS sampling sites in Lake Erie, Green Bay, etc. 
which were associated with the collection of large volume samples for 
the NASA/LaRC optical parameter analysis. Results from these laboratory 
analyses are reported in Appendix A. 
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4.0 


GEOMETRIC CORRECTION METHODS 

During phase I of this projejct a scanner model was developed which 
permitted the transformation of CZCS line and point coordinates to earth 
latitude and longitude. CZCS scanning geometry includes variable tilt 
angle for Great Lakes viewing combined with ground control points to 
derive the appropriate transformation matrix. Existing mapping software 
was modified to accommodate the CZCS frame size and mapping polynomials. 
Geometric correction of an image results from two operations. First 
mapping polynomials are generated to define the transformation from raw 
original image to the corrected image. Second, the corrected image is 
created from the uncorrected image using these mapping polynomials. Two 
fifth-order, twenty-one term polynomials were used in this process, one 
for each dimension of the image. These polynomials defined the trans- 
formation which makes the corrected image conform to a given map pro- 
jection as well as adjust for viewing distortions such as satellite 
position, satellite motion, and earth motion. 

The scanner model developed for CZCS is based upon the ERIM experi- 
ence with Landsat and in its present form takes each image control point 
in turn and proejcts it to the earth's surface. The line number for 
each point is used to interpolate for the latitude, longitude, and 
altitude from the values supplied with the tape reference data for each 
scene. The point number is used to calculate the mirror scan angle 
which together with the reported tilt angle determines the scanner line 
of sight vector in spacecraft coordinates. A series of single angle 
rotations through the angles of roll, pitch, heading provide the 
transformation to earth centered coordinates. 

There are two basic uses for the derived scanner model. First, the 
image control points and their corresponding map control points may be 
easily evaluated for consistency with other points and any outliers 
rejected. Second, the coefficients for global mapping polynomials used 
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in the resampling process can be derived by a fit to the model rather 
than to the points themselves. 

A number of scanner variables including satellite altitude and 
longitude are converted to line and pixel location in the resulting 
image via the equations which define the desired map projection. The 
process used to resample a corrected image works in the reverse direc- 
tion. For each pixel location in the resulting image the software 
calculates the corresponding location in the original image. Operation 
of the scanner model in reverse so as to select original pixel locations 
which correspond to a given location in the corrected image is most 
difficult. Alternatively the pair of twenty one term polynomials is 
generated to satisfy this mapping requirement. One polynomial describes 
the east-west position and the other the north-south position. The 
selected set of image control points are used to generate coefficients 
for each term in the polynomial. The complete set of coefficients 
defines the mapping from the original image to the correct projected 
image. 

4.1 RESAMPLNG OF THE CZCS 1980 IMAGERY 

Each of the five scenes collected during the summer 1980 CZCS Vali- 
dation Experiment as listed in Table 3.1, were resampled with the same 
polyconic projection space with the central meridian set of 84.0 
degrees. Thus, each resampled image was coregistered to within the scene 
to scene resampling errors. A data file with one hundred and twenty- 
four control points was compiled for the Great Lakes region. Most of 
these points could be easily identified in the CZCS imagery and were 
somewhat evenly distributed throughout the region. Control points were 
selected from this data file for each scene which could be identified in 
the raw image. In the case of the June 11 and July 23 scenes (orbits 
8243 and 8823) most of the Great Lakes were visible and multiple frames 
(two or three) from the same orbit were required to cover the entire 
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region (see Figure 4.1). In this case resampling polynomials were gen- 
erated for the frame with greatest coverage and extended to apply to the 
entire data set. The resampling software did not permit generation of 
mapping polynomials based upon a variable size input space. The possi- 
bility existed to independently generate mapping polynomials for each 
frame but the resulting problems of boundary matching made the approach 
impractical. Results from independent polynomial mapping compared well 
with that obtained by domain extension. 

All of the CZCS frames were resampled with a 500 meter pixel size, 
but the pixel size in the input space, of course, varies with scan 
position. Radiometric resampling employed a nearest neighbor algorithm. 
Spatial accuracy of this resampling procedure is thus a function of scan 
position and ground point control. The approximate position of the 
nadir track in the resampled space for each of the scene date is shown 
in Figure 4.1. In some cases only a few control points could be applied 
since most of the Lakes were cloud covered. The lack of control points 
was generally compensated by selecting a smaller area for resampling and 
one where the resampling polynomial was considered to be representative. 
The worst case occurred with the July 17 date where only a portion of 
Lake Erie containing the surface sampling site was cloud free. The rest 
of the Great Lakes region were essentially cloud covered. This condi- 
tion coupled with the large scan angle produced the largest resampling 
RMS errors (approximately 5000 meters E-W and 2100 meters N-S). Typical 
resampling errors were estimated to be 200 meters in the N-S direction 
and 1875 meters in the E-W direction [3]. Since the scanning direction 
of CZCS is nearly east and west at this latitude the RMS east-west 
errors were found to be much larger than the north-south errors. 

Subsequent radiometric analyses of the resampled CZCS data sets 
were made on individual pixels corresponding to ship sampling sites and 
on subscenes containing a single Great Lake or sublake (e.g.. Green Bay, 
western Lake Superior, etc.). The individual pixels were extracted 
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based upon the lat/long line/point relationship implied by the polyconic 
projection space and the reported ship position at the sampling site. 
Because of the resampling errors and the high variability of lake 
waters, some of these pixels are more representative of the radiometry 
of the sampling position then others. 
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DEVELOPMENT OF AN ATMOSPHERIC CORRECTION ALGORITHM 
FOR THE GREAT LAKES 

Before CZCS data can be effectively used to distinguish between and 
monitor coastal water types, methods are needed for removal of spatial 
variations in aerosol path radiance which dominate radiance measurements 
made by the satellite sensor. In the present study efforts were focussed 
on the comparison and evaluation of existing atmospheric algorithms for 
use with Great Lakes data and the development of improved algorithms. 
Algorithms were sought which would have minimum requirements for surface 
truth and which would produce scene independent results. The abilities 
of three existing algorithms and one developed under the present study 
to quantitatively remove haze effects when applied to Great Lakes scene 
data where evaluated and compared. 

The atmospheric algorithms considered are (1) the Gordon Algorithm 
[4], (2) the Smith and Wilson iterative algorithm [5], (3) the pseudo 
optical depth method [6], and (4) the residual component algorithm. 
Each of these algorithms was adapted to extract required scene dependent 
parameters during an initial pass through the data set. The performance 
of each algorithm was evaluated on the basis of four comparisons: (1) 
derived upwelling radiance spectra with water radiance model calcula- 
tions, (2) derived aerosol path radiance spectra with atmospheric model 
predictions, and (3) predicted radiance in the form of color coded 
images with uncorrected data. These latter maps show the extent that 
spatial features of the aerosol haze were removed by each algorithm. 

Atmospheric path radiance is produced by a combination of Rayleigh 
and aerosol scattering. The amount of path radiance generated by 
Rayleigh molecular scattering is a function of viewing and solar geome- 
tries and can be calculated directly from an analytical model. The 
radiance due aerosol particles varies spatially and cannot be predicted 
by analytical means unless local atmospheric optical thicknesses are 
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well known. Because of spatial variations in the concentration of aero- 
sol particles it is necessary to use an atmospheric correction procedure 
that applies corrections on a pixel-by-pixel basis. Effective correc- 
tion algorithms of this type have been demonstrated by Gordon, et al [7] 
for the open ocean. The ocean waters for which these algorithms work 
are classified as type I by Morel and Prieur [9]. These algorithms will 
not work in the more turbid case II coastal waters. The waters of the 
Great Lakes are by in large of type II but cover a wide range of optical 
properties and concentrations of suspended particles. Waters of central 
Lakes Superior and Huron approach those designated as type I. 

5.1 RADIATIVE TRANSFER PROCESS IN THE LAKE ATMOSPHERE SYSTEM 

The components of the radiative transfer model for a lake atmos- 
phere system are described by the equation below. 

L( X) = lJ(x) + LrU) + L a (X) + L S (X) + L U (X) . T( X ) (2) 

where L(x) is the radiance measured at the satellite; L*(x) is the 
Rayleigh path radiance as generated by molecular backscattering of solar 
irradiance; L^(x) is the Rayleigh path radiance due to scattering of 
irradiance after it has reflected from the surface toward the sensor, 
L a (x) is the aerosol path radiance, L^(x) is the reflected sky radiance, 
Ly(x) is the upwelling radiance from beneath the water surface, and T(x) 
is the atmospheric transmittance. 

The Ly(x) term contains the information on the water constituents 
and may only represent a small fraction of the total signal received at 
the satellite. The Rayleigh scattering phase function and the sensor 
view geometry can be used together to directly calculate the Rayleigh 
components at each pixel. Unlike the Rayleigh scattering process it is 
not possible to approximate the aerosol path radiance explicitly since 
aerosol concentration is highly variable. In principle, however, the 
radiance added by aerosol scattering could be calculated if the concen- 
trations and optical properties were known throughout the scene. Each 


36 


2p 

algorithm must, therefore, use the data itself to derive the aerosol 
component at each pixel in order to subsequently estimate the Ly term. 
This approach relies on making two assumptions. The first of these 
assumptions is that the band ratio of aerosol path radiances is 
independent of pixel position within the scene. This statement is 
equivalent to saying the aerosol concentration varies from 
pixel -to-pixel but aerosol type remains constant. This assumption is 
largely supported by the fact that the aerosol scattering phase function 
is essentially independent of wavelength. In addition, it is necessary 
to assume some characteristic of the upwelling radiance spectra which is 
invariant over the entire scene. 

The approach taken to calibrate the algorithm involves collecting 
sea truth data at one position within the scene. First, the downwelling 
irradiance or atmospheric optical thickness is measured at each of the 
CZCS wavelengths. Second, the upwelling radiance just below the surface 
is measured directly or model estimated from the measured chlorophyll 
pigment concentration. Together these two measurements can be used to 
estimate the aerosol path radiance ratios. The atmospheric correction 
is then adjusted over the calibration region so that the satellite 
derived upwelling radiances (or reflectances) match the surface measured 
or model estimated values. Gordon has refined this approach for studies 
of the open ocean with his concept of normalized clear water radiances 
[15]. For this study normalized water radiances for Great Lakes waters 
were developed for a Lake Superior calibration area. 

5.2 SOLAR IRRADIANCE AND ATMOSPHERIC TRANSMITTANCE 

The mean extraterrestrial solar spectral irradiance, E Q , has been 
determined experimentally for each Julian Day. Values of E 0 for each of 
the CZCS wavelengths are given in Table 5.1. The solar irradiance is 
attenuated during the passage through the atmosphere caused by interac- 
tion with air molecules and aerosols during which photons are either 
scattered or absorbed. The quantity of radiance transmitted depends on 
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TABLE 5.1 ATMOSPHERIC OPTICAL PARAMETERS 


Wavelength (nm) 


443 520 550 

Extraterrestrial Irradiance (mW/cm . sr) 

186.42 185.34 184.76 

Rayleigh Optical Thickness 

0.2316 0.1224 0.0964 

Ozone Optical Thickness 

0.0069 0.0237 0.0390 

Initial Aerosol Optical Thickness 

0.320 0.300 0.280 


670 

151.52 

0.0442 

0.0226 

0.250 


the optical thickness and the slant path. The optical thickness charac- 
terizes the attenuation properties in a vertical path and can be 
separated into molecular and aerosol optical thicknesses as given in the 
equation below. 


x(x) - Xp(X) + X^(X) + T Q3 ( -^ ) 


(3) 


where x (x) Q3 is the optical thickness of the ozone layer, x R (x) is the 
Rayleigh optical thickness due to molecular scattering and x^(x) is the 
aerosol optical thickness due to scattering and absorption by aerosols. 
The ozone optical depth varies with season and latitude and the values 
used for this study are given in Table 5.1. Because the quantity of 
molecules in the atmosphere is essentially constant, the Rayleigh 
optical thickness varies approximately inversely with the fourth power 
of wavelength. The aerosol optical thickness, on the other hand, varies 
throughout the scene. 
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The diffuse transmittance, Tp(x), is usually applied to components 
of radiance between the water surface and the sensor. This 
transmittance is given by 


T D (X) = EXP[-(a R . x R (x) + t 03 (x) + a A . t a (x))/m] (4) 

where y is the cosine of the viewing zenith angle. The diffuse 
transmittance is essentially the same as the beam transmittance except 
that each component of the optical thickness has been multiplied by a 
factor to account for scattering of radiance from adjacent pixels which 
are scattered into the sensor field of view. The ozone layer absorbs 
and does not scatter photons so there is no change for the diffuse 
calculation. The Rayleigh scattering phase function p(ijj) is symmetrical 
as given below and, therefore, one half of the photons scattered from 
the beam will be replaced by those scattered into the beam and a factor 
of one half is, therefore, appropriate to this calculation. 

p(ip ) = 3/4(1 + cos 2 (ip ) ) (5) 

The aerosol scattering is asymmetric with a large forward component 
and, therefore, the diffuse transmittance is only weakly dependent on 
the aerosol optical thickness. A coefficient of 0.1 has been proposed 
but 0.0 can be effectively used where haze levels are very low. 

5.3 CALCULATION OF THE RAYLEIGH PATH RADIANCE COMPONENTS 

The radiative transfer equation can be solved using the Rayleigh 
scattering phase function to derive closed form solutions for each of 
the components. Only single scattering is considered since multiple 
scattering is negligible for this case. The Rayleigh path radiance for 
directly scattered radiance is given by the following equation. 
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y + m 0 


1 E o( x ) 

L^IXsPjUq) — • [1 ~ 6Xp(- T ^ . )] . Tq^X, y, Up) 

4 * (y + Uq ) y • y Q 

where T Q3 (x, y, y ) is the transmittance for a combined downward and 
upward trips through the ozone layer, y and y Q are cosines of the view 
and solar zenith angles and Eg is the extraterrestrial irrradiance. The 
scattering angle (as shown in Figure 5.1) is given by the following 
cosine equation where 9 and 9 o are, respectively, the view and solar 
zenith angles and 0 and 0 Q are the view azimuth and solar azimuth angles 
to the solar plane. 


( 6 ) 


cos(^_) = -cos 9 . cos 9 q + sin 9 . sin 9 Q . cos (0 - 0 ) (7) 

This approximate solution is the quasi -single scattering solution which 
is discussed in section 5.5. It has been assumed by some investigators 
that multiple scattering is negligible in which case the above solution 
can be approximated by 

E 0 (X) P(*J 

L^( X ,y ,y Q ) - . Tg . Tgg( X ,y ,y Q ) (8) 

4n y 

which is essentially the double delta solution [17]. 

2 

The Rayleigh component Lj^(x) of direct radiance which is first 
reflected from the water surface and then scattered to the sensor is the 
same as the equation above, but includes a Fresnel factor for the solar 
zenith angle. The scattering angle is given by the following cosine 
equation. 

cos(i^ + ) = cos 9 . cos 9 q + sin 9 . sin 9 Q . cos (0 - 0 Q ) (9) 
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FIGURE 5.1 GEOMETRY OF RAYLEIGH SCATTERING ANGLE 



FIGURE 5.2 CZCS VIEWING GEOMETRY IN THE RESAMPLED IMAGE 
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The last term accounts for the photons which are scattered on the 
downward side and reflected directly to the sensor. The equation for 
this term is the same as above, but where the solar Fresnel factor is 
replaced with one for the view angle. These three components together 
constitute the Rayleigh path radiance, and were included in each of the 
algorithms tested. 

In making these calculations for resampled CZCS data one must 
relate the resampled pixel coordinates to the sensor viewing geometry. 
This relation is most easily acquired by locating at least two nadir 
points in the resampled scene. The nadir track can be assumed to be 
approximately linear with line number in the resampled space (see Figure 
5.2). Because of rotation, the track will have angle z, with the 
geographically orientated resampled space. For a sensor tilt angle X, 
the view angle 9, and view azimuth i are given by the equations below. 

9 = tan -1 [tan 2 x + (aP + aL* tan c) cos c/H 2 ] 1//2 


<t> = tan _1 [(AP + AL-tan C) cos C/H-tan x] - 5 - 0. 


( 10 ) 


where S is the angle between the nadir track, and the resampled scan 
line, aP and aL are the point and line offsets (in km) from a reference 
point on the nadir track, and H is the spacecraft altitude. 

5.4 CALCULATION OF AEROSOL PATH RADIANCE 

Theoretically we can calculate the aerosol path radiance components 
in the same as those for Rayleigh scattering if we knew the aerosol 
optical thickness, absorption coefficient (or single scattering albedo) 
and the scattering phase function. Since the aerosol optical properties 
are highly variable in both time and space such analytical calculations 
are not possible for algorithm calculations. Thus the calculation must 
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be accomplished by some indirect means which will be the focal point of 
each of the correction schemes discussed below. 

5.5 ATMOSPHERIC WATER RADIANCE MODEL 

The radiance at any point in the lake atmosphere system can be 
obtained in principle by the solution of the radiative transfer equation 
with proper knowledge of critical optical properties. A full descrip- 
tion of the optical properties of a medium include the absorption 
coefficient and volume scattering function. An exact solution is not 
obtainable in close form for this integral equation even for the 
simplest geometry. Under this circumstance two approaches are possible. 
First, we can develop an exact numerical solution with Monte Carlo simu- 
lation techniques. This method is especially powerful and can incorpor- 
ate spatially varying optical properties. The chief disadvantage to 
this approach is the large number of iterations required to obtain a 
solution for a single case. Further, a large number of cases need to be 
examined in order to gain proper insight into the problem. A more 
practical approach and the one taken here, is to develop approximate 
solutions which give good results over a range of conditions. Insight 
can sometimes be gained by merely examining the functional form of the 
approximate solution. 

Using this approach radiative transfer model was developed to 
predict the CZCS radiances at the top of the atmosphere which includes 
photon interaction in the atmosphere, through the water surface and 
beneath the water surface. A flow diagram for this model is shown as 
Figure 5.3. The model was designed to accept specific view and solar 
geometries as well as aerosol type characteristics giving it the 
capability to simulate any selected CZCS pixel radiance. 

The atmospheric portion is based upon a previous model developed by 
Turner et al [21]. The extraterrestrial irradiances, Rayleigh optical 
thicknesses, and ozone optical thicknesses are the same as reported 
above. The Turner code was modified to include the ozone layer, the 
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FIGURE 5.3 WATER RADIANCE MODEL FLOW CHART 
























reflected diffuse sky radiance, and the reflected direct irradiance 
which is scattered to the sensor. The aerosol phase function was based 
upon an Elterman type atmosphere [22]. Turner used the following 
aerosol scattering phase function which allows direct solution of the 
radiative transfer equation. 

p(u, 0; VI 1 , 0') = 4irF . y(y - y') y(0 “ 0‘ ) + 4irB . y(y + y') y(0 - 0' + it) (11) 

where F and B are the forward and backward scattering constants and y is a 
delta function. The solution was re-entered into the transfer equation to 
get a second order solution in what is referred to as the double delta 
method. The Turner solution has been compared to exact Monte Carlo 
solutions by Lyzenga [23] and found to give acceptable results in the 
forward direction but overestimated the radiance in the anti -solar direction 
and underestimated the radiance in other directions. The double delta 
solution thus gives a good solution for the downward irradiance and the 
reflected upward scattered direct radiance. The double delta solution will, 
in general, work poorly for strongly forward scattering functions which are 
characteristic of aerosols, but give good results for di symmetric Rayleigh 
scattering. These solutions were then incorporated into the present model. 

A quasi-single scattering solution was used to estimate the radiance that is 
backscattered to the sensor during the downward leg. The quasi-single 
scattering solution uses the following phase function which includes the 
actual aerosol scattering function in the backward direction as given by the 
following equation. 

P ( y , 0; u'» 0') = 4nF . Y(y-y‘ ) Y ( 0-0 1 ) + p(y, 0;y , ,0 , )yy l <.0 (.12) 

The quasi -single scattering solution yields a hemispherical reflectance 
which agrees well with exact calculations. While this solution estimates 
the backward scattered radiance it does not predict forward scattered 
radiance. In addition, the errors in the backscattered radiance increase as 



the single scattering albedo approaches one (i.e., a highly turbid 
atmosphere). The quasi-single scattering solution does account in part for 
multiple scattering of photons. The Turner model was also modified to 
include the diffuse transmittance as defined previously for the component 
upwelling from beneath the surface and for the reflected diffuse sky 
radiance. The diffuse transmittance also partially accounts for possible 
multiple scattering effects. 

In developing a water radiance model, it is generally assumed that 
the inherent optical properties arer linear functions of concentration 
of each constituent in the medium. Under these circumstances the total 
absorption, backscatter, and total scattering coefficients can be 
expressed by the following relationships. 

N 

a = a + z a -c . 
w i =1 * 1 


Bb = Bb + z Bb.c. ( 13 ) 

w .,11 
1=1 

N 

b ■ b W + Vi 

where a w , Bb w , and b w are the absorption backscattering and total scat- 
tering coefficients for pure water and a^ , Bb^ , and b^ are the absorp- 
tion, backscattering, and total scattering cross sections, respectively 
for constituent i with concentration c^. These relationships are ap- 
proximations since the actual optical properties depend on other factors 
such as particle size distribution. Having made these approximations a 
useful relationship can be made between these inherent optical 
properties and the apparent optical properties of the medium including 
volume reflectance. The water reflectance model used for this study was 
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that developed by Gordon, et al . [8] which is commmonly referred to as 
the power series approximation and is given below. 


R~(X) = 0.0001 + 0.3244X + .1425X 2 + .1308X 3 (14) 

where x is Bb/(a+Bb) and R~ is the subsurface hemispherical reflectance 
and the (-) signifies a subsurface component. The coefficients were 
obtained by fitting the power series to Monte Carlo calculations of the 
subsurface water reflectance for solar zenith angles less than thirty 
degrees. These calculations were made for hemispherical irradiance 
reflectance and for this model we assume a uniform reflectance distribu- 
tion so that the subsurface radiance reflectance can be estimated. 
There have been some studies to indicate that the relationship between 
these quantities is somewhat greater than it [8]. However, for this 
study it was assumed that the radiance reflectance r can be approximated 
as 


r‘(u, u Q ) = R^(u 0 )/* ( 15 ) 

The above surface upwelling radiance from beneath the surface can then 
be calculated using the following relationship. 


r + (y, y Q )= — r ”(u, m q ) ( 16 ) 

2rm 

where T^ and T£ are the surface transmittances ( i . e , . 1 - p, p = Fresnel 
reflectance) for the incoming and outgoing light. 

Derivation of optical properties for the water reflectance models 
are discussed in section 6.1. 
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5.6 EXISTING ATMOSPHERIC ALGORITHMS 

Three existing algorithms were selected for comparison and evalu- 
ation with CZCS data collected over the Great Lakes. These included 
(1) the NOAA NET algorithm developed by Gordon et al . [4], (2) the 
iterative algorithm of Smith and Wilson [5], and (3) the pseudo-optical 
depth algorithm developed for the Great Lakes by Jain et al . [6]. Each 
of these algorithms was adapted to operate in two modes. First, the 
algorithm can be used in the calibration mode to extract necessary 
aerosol parameters from a hazy portion of the scene. Second, the 
algorithm can, with these data derived parameters, process the scene on 
a pixel -by-pixel basis for removal of aerosol and Rayleigh path 
radiance. 

1. Gordon Algorithm 

The NOAA atmospheric correction algorithm for CZCS developed by 
Gordon assumes the upwelling radiance at 670 nm is zero for all pixels 
to be corrected. This assumption has been shown to be valid for ocean 
waters and using the aerosol ratio values the direct calculation of t*L w 
by way of the following equation. Here L w is the sum of the upwelling 

W 

radiance L u plus the reflected sky radiance. 

L(A.) = L r (\.) + T D (x i )*L u (x i ) + k ij .(L(X i ) - L R (x i ) ) (17) 

where k^- = e (x 1 -, x J )[F 0 (x 1 )/F 0 (x J )] and F Q is the instantaneous extra- 
terrestrial solar irradiance reduced by the two way transmittance 
through the ozone, and c is the ratio of . oj and the aerosol 

scattering phase function at the two wavelengths, x i and Xj. In this 
manner one has separated out the factor e which does not change with 
respect to pixel position. Once the e value is determined for a 
particular scene the aerosol path radiance can in principle be estimated 
for each pixel and subsequently the water leaving radiance. The 
determination of £ is accomplished by having at least one pixel in the 
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scene where the value of L has been measured or estimated using the 
concept of clear water radiances [24]. In equation 17 k. . is the band 

’ J 

ratio of the aerosol path radiances. 


2. Smith and Wilson Algorithm 

In coastal waters the subsurface upwelling radiance at 670 nm is 
non-zero largely because sediment concentrations are relatively high. 
Smith and Wilson [5] developed an iterative procedure for these waters 
which is based upon known spectral characteristics and which allows L 

w 

at 670 nm to be non-zero. The algorithm utilizes the following 
relationship determined from measured upwelling radiance data. 

L w (670 nm) = 0.0829 l_ w (443)-R l3 -1 * 661 (.18) 


where R 13 = L w (443 nm)/l_ w (550 nm). Initially one assumes that L w = 0 
and uses the above equation to calculate the upwelling radiance at each 
of the other wavelengths per the NET algorithm. The values of L w at 443 
and 550 nm are then substituted back into the above equation and a new 
estimate of L is obtained. This process is repeated until there is 
acceptable convergence to constant values of L.,(a). 

W 

3. Pseudo-optical Depth Algorithm 

In this algorithm path radiance is assumed to be a linear function 
of surface albedo for a specific aerosol optical depth. The discrete 
ordinate method (DOM) solution to the radiative transfer equation is 
used to derive the linear relationships between surface albedo R, and 
total path radiance Lp including the aerosol and Rayleigh components. 


L p = C 1 + C 2 ’ R w 


( 19 ) 


The computation assumes the ozone and Rayleigh standard atmosphere of 
Elterman and a Deirmendjian haze (L) continental aerosol [21]. The 
Elterman aerosol optical depth is scaled according to the experimentally 
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derived values. Absorption was assumed to be small and not included in 

these calculations. For fixed values of the surface albedo and view 

angle these calculations were fitted to a second order polynomial in the 

aerosol optical thickness. Accuracy of the computed Lp with the 

approximate equations was reported by Jain [6] to be within 1.1% of the 

actual DOM calculations. A linear equation was used to describe the 

relationship between L and the solar zenith angle for R < 0.25 and 

p w — 

with cos 9 for 0.25 < R w < 0.50. These relationships together provided 
a means to determine L p from R w , 9 Q , and x a * In addition, these 
relationships can be inverted to approximate from a measured L p , R w , 
and 0 Q . The algorithm procedure is as follows. For each pixel the path 
radiance is estimated at the reference wavelength (670 nm or 750 nm) and 
the upwelling radiance is assumed to be known. A pseudo aerosol optical 
depth is then estimated from the DOM based relationships. 

The ratio between this pseudo optical depth to a previously 
measured optical thickness at some location within the scene is used to 
estimate by proportion the aerosol optical depth at each of the other 
wavelengths. The underlying assumption here is, of course, equivalent 
to that common to all algorithms as discussed above. These aerosol 
optical depths are then used with the derived relationships to calculate 
the path radiance as a linear function of the surface albedo. The 
surface albedo can be estimated as 


* • L W U) 

R w (A) = E q (x) - T 0 (X) 
In this case the satellite radiance is 


L(t) - L w (x) . T d U) + Cj + • it • L w (x) 


( 20 ) 


(. 21 ). 
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where Tq is the downwelling beam transmittance which can be calculated 

using the previously estimated aerosol optical depths. The above 

equation can now be solved for L which does include the radiance from 

w 

reflected diffuse sky irradiance. 

5.7 DEVELOPMENT OF A NEW CORRECTION ALGORITHM FOR THE GREAT LAKES 

A fourth algorithm considered in this study was centered around the 
need to extract the needed algorithm parameters from the data itself. In 
this approach the characteristic CZCS spectra are used to develop a 
diagnostic for haze type which is directly related to the previous 
approach but which avoids the necessity of assuming zero upwelling 
radiance in the 670 nm band. Basically the algorithm assumes that over a 
uniform small patch of water the upwelling radiance will be essentially 
constant as will the Rayleigh components of path radiance. In this case 
the pixel -to-pixel variations are due solely to the variations in the 
aerosol path radiance components. The orientation of this radiance 
vector can be determined from the scene data independent of actual 
measurements of the upwelling subsurface radiance. The aerosol vector 
direction cosines are obtained as the eigenvector of the first principal 
component. These cosines are related directly to the k-ratios of the 
aerosol path radiances by 


cos Tj = (1 + z k^. 2 )" 1 ^ 2 ( 22 ) 

j=l 

j^i 

where Y . is the angle between the aerosol vector and the aerosol path 
radiance axis for band i. The orientation of the aerosol vector is 
assumed to be constant throughout the scene or portion to be atmospheri- 
cally corrected. This assumption is equivalent to assuming a constant 
aerosol type with the length of the aerosol vector allowed to vary from 
pixel -to-pi xel . Example first component eigenvectors are listed in 
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Table 5.2 for the Great Lakes scenes with the percent of total variance. 
Our analysis of these results suggest that a good estimate accounts for 
over ninety percent of the variance. In cases where the aerosol path 
radiance components are small the eigenvector is not uniquely determined 
and surface measurements of the optical thickness and/or measurements of 
the upwelling radiance are required to estimate proper eigenvector 
components. Like algorithms one and two the residual component algor- 
ithm relies on a scene invariant relationship between the upwelling 
radiance in each band such that 

F(L u (x i ). 1 = 1 ,2,3,4) = 0 ( 23 ) 

In principle we will not know the length of the aerosol vector at 
any pixel and, therefore, must estimate this value with the algorithm. 
The water radiance vector, hopefully is always oriented in some other 
direction from that of the aerosol vector. Previous model simulation 
studies have indicated that 18.7 to 26.1 degrees separate the aerosol 
vector and suspended minerals and chlorophyl 1-a vectors respectively 
[3]. There is, in any case, a large component of the water radiance 
vector cooriented with the aerosol vector. A relationship was, 
therefore, developed between the projected modulus (brightness) of the 
water . radiance vector, | L (x) |, and a spectral characteristic R* which 
is invariant throughout the scene. The spectral characteristic selected 
for the residual component algorithm was the following combination of 
upwelling radiances. 


R*(L W U)) = 


L w (550 nm) . L w (520 nm) 


L w (443 nm) 


(24) 


Figure 5.4 shows the water radiance model derived relationship for a 
June 11, 1980 aerosol vector orientation. The water radiance values were 
calculated using the Great Lakes optical model described in section 6.3. 
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TABLE 5.2. FIRST PRINCIPAL COMPONENTS DERIVED FROM CZCS DATA SETS 
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R* 


Figure 5.4 R* vs. brightness (modulus) of the upwelling 

radiance spectra using simulated cases generated 
with the Great Lakes optical model. 
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For small values of R* the L (550 nm) . L ( 520 nm) product ratio is 
correspondingly small and the water radiance brightness is dominated by 
the contribution from the 443 nm band. Waters in this case have very 
low concentrations of chlorophyll and are essentially free of suspended 
minerals and dissolved organics. As R* increases so do the constituent 
concentrations. Surprisingly the corresponding water radiance bright- 
ness values decreases, passes through a minimum, and then increases 
asymmtotally to a maximum. These data imply a higher order polynomial 
in R* would describe the relationship. A study by Austin and Petzold 
[25] of a large number of ocean upwelling radiance measurements showed 
that the radiance in each band can be expressed as a power function in 
R* suggesting as in the present study that some higher order polynomial 
would describe the vector modulus of the water radiance. For the 
residual component algorithm the following linear approximation was used 
to describe this relationship. 

I L (X) I = 0.517 - 0.097 R*( L ( X ) ) R* < 1.4 
w w 

(25) 

I Ljx) I = 1.318 + 1.17 R*(L (x) ) R* >_ 1.4 

w w 


In applying the residual component algorithm the Rayleigh path 
radiance components are first calculated and subtracted from the total 
satellite measured radiance. The Rayleigh components are estimated here 
using the double delta solution rather than the QSS approximation. Next, 
the aerosol path radiance is estimated at 443 nm and the previously 
derived a-ratios are used to estimate the corresponding path radiance at 
each of the other wavelengths. These radiances are then subtracted at 
each wavelength from the balance and the length of the residual compon- 
ent vector as projected on the aerosol vector is calculated. The 
residual radiances are used to also calculate R* and make a model 
predicted estimate of the brightness value. The two moduli are then 
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compared and the length of the aerosol vector is modified so as to 
reduce the difference in these two estimates. In effect we are adding 
or subtracting the amount of aerosol haze present while keeping the type 
constant. The procedure is conceptually illustrated in Figure 5.5. Two 
or three iterations are usually sufficient for convergence to a constant 
set of residual component radiances which are then equal to 

T d U) • L W (X). (26) 

5.8 CONVERSION OF DIGITAL COUNTS TO RADIANCE VALUES 

Since each of the above algorithms operates on CZCS measured 
radiances the computer code for each algorithm contained a section to 
convert count data to radiance values. These conversion relationships 
are critical to obtaining any meaningful results where a quantitative 
approach has been made to understand the radiometric components. 

The necessary radiometric calibration constants are nomially re- 
ported in the header of each NOAA processed CCT. However, NET calibra- 
tion studies have shown that use of such constants would produce incor- 
rect radiance values primarily because these constants do not account 
for sensor deterioration. 

Two independent investigations were made by H. Gordon and R. Austin 
et al [26] into the calibration problem using repeated CZCS measurements 
over predictable very clear waters such as those of the Sargasso Sea. 
In both studies the objective was to start with the original prelaunch 
gain correction constants and determine an adjustment factor G(x,j) for 
each orbit k. The sensor radiances L ( A ) are then related to the digital 
counts N(A,j) by 


L(A) = (A(a) . N(A,j)K) + B( A) )G( A, j ) (27) 
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where A and B are the gain correction constants. Unfortunately, the two 
investigations, which used tens of scenes over deep ocean targets with 
uniform upwelling radiance, came up with different results. The two 
studies produced two fits of G with orbit number with different mathe- 
matical form. The G values were found to be similar at some orbit 
numbers and quite different at others. The G calibration factors used 
in this study were based upon a set of formulas developed by H. Gordon 
[27] which are listed with the gain 1 and 2 correction constants in 
Table 5.3. Also given are the G values for orbits 8243, 8699, 8740, 


8823, and 9017 (June 
respectively. 

11, July 14, 

July 17, 

July 23, 

and August 


TABLE 5.3 

RADIOMETRIC 

CALIBRATION CONSTANTS 

1980 Dates 

Orbit(j ) 

443 nm 

G(\,j) values 
wavelenth (a) 

520 nm 550 nm 

620 nm 

June 11 

8243 

1.189 

1.022 

0.969 

1.00 

July 14 

8699 

1.201 

1.025 

0.969 

1.00 

July 17 

8740 

1.202 

1.025 

0.969 

1.00 

July 23 

8823 

1.204 

1.025 

0.969 

1.00 

August 6 

9017 

1.209 

1.027 

0.971 

1.00 

A( X ) 

(Gain 1) 

.0445 

.0310 

.0247 

.0114 

B(X) 

(Gain 1) 

.0396 

.0636 

.0799 

.0114 

A(x) 

(Gain 2) 

.0359 

.0249 

.0202 

.0897 

B ( X ) 

(Gain 2) 

.0528 

.0883 

.0625 

.0359 


Above Calibration Constants from H. Gordon [26]. 

5.9 ALGORITHM CALIBRATION 

Each of the existing atmospheric correction algorithms and the 
residual component algorithm developed as part of the present study were 
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applied to the June 11, 1980 scene for purposes of comparison and 
evaluation. Large quantities of aerosol haze were present over Lakes 
Michigan and Superior on this date. These aerosol concentrations were 
highly variable and easily identifiable in each CZCS band especially 
over a major portion of Lake Superior where conditions ranged from clear 
to haze sufficiently heavy to optically obscure the upwelling radiance. 

The algorithms were implemented by making two passes through the 
data set. In the first pass the necessary scene dependent parameters 
were estimated for each algorithm as discussed below. The second pass 
removed the aerosol components of path radiance on a pixel -by-pixel 
basis and transformed the corrected radiances to a nadir viewing 
geometry. The results of this algorithm processing were summarized in 
products which consisted on numerical transects and color coded radiance 
maps. 

The four algorithms were individually calibrated for the Lake 
Superior subscene using a block of 100 pixels located north of 
Michigan's Keewanaw Peninsula in the central part of the Lake (see 
Figure 3.4). Typical chi orophyl 1 -a pigment levels have been reported by 

3 

Munawar et al . [28], 1978 to average 1.1 mg/m for all but the nearshore 
waters with an expected lower seasonal value for the calibration area of 

3 

0.8 mg/m . Chlorophyll levels are very stable for this part of Lake 
Superior and perhaps the most predictable for all of the lakes. 
Year-to-year productivity cycles in Lake Superior are not expected to 
change dramatically based upon year-to-year measurements made at the 
same time over a period of several years, including 1980, in the open 
waters of Lake Huron [20]. While no surface truth data were collected 

3 

from this area on June 11, 1980 the estimated concentration of 0.8 mg/m 
is considered to be within chlorophyll-a measurement error of actual 
values. Surface truth collected on this date were collected in Lake 
Michigan near Grand Haven, Michigan. This area displayed little 
radiometric effects from haze and had highly variable concentrations of 
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chi orophyl 1 -a and suspended mineral particles associated with the 
discharge from the Grand River. For these reasons the Grand Haven 
surface truth data could not be used to calibrate the algorithms to be 
applied over the whole of Lakes Superior and Michigan. However, the 
optical thickness measurements were used to estimate the a-ratio 
coefficients which were considered to be valid for much of the southern 
portion of Lake Michigan. 

Estimation of the upwelling radiance for the Lake Superior calibra- 
tion area was accomplished using the water radiance model discussed 
above and in section 5.5. A normalized upwelling radiance was estimated 
for Lake Superior which was independent of solar elevation and view geo- 
metry. These radiance values were then used to estimate the upwelling 
radiance from beneath the surface at the calibration area. The viewing 
geometry and normalized radiances are shown in Table 5.4. In the case 
of the Gordon algorithm the average satellite radiance, estimated 
upwelling water radiance, and average calculated Rayleigh path radiance 
were used to estimate the average aerosol path radiance at each 
wavelength and subsequently the epsilon constants which are related to 
the k-ratios and which account for scene-to-scene differences in aerosol 
type. A similar set of epsilons were derived for the Smith and Wilson 
iterative algorithm but were determined to be slightly different since 
the upwelling radiance at 670 nm was not assumed to be zero. In the 
case of the pseudo optical depth algorithm the average path radiance 
over the calibration area was used to estimate the corresponding aerosol 
optical depth from the DOM derived relationships. These aerosol optical 
depths were then used as described above to estimate the aerosol optical 
depths at other pixels within the scene. For the residual component 
algorithm the calibration area data were used to derive a first 
principal component eigenvector containing the direction cosines of the 
aerosol vector. The model generated upwelling radiance values were not 
needed for this latter algorithm calibration. The scene dependent 
parameters derived for the four algorithms are given in Table 5.5. 
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TABLE 5.4 VIEWING GEOMETRY FOR LAKE SUPERIOR CALIBRATION AREA 
Solar zenith angle: 25.4 degrees 

Solar azimuth: 152.4 degrees 

Average view angle: 25.6 degrees 

Sensor tilt angle: 20.0 degrees 

Calibration area size: 100 pixels 


Normalized radiances L (mW/cm.sr.y) (1. . .4) = (1.495, 0.980, 0.725, 
0.042) 


TABLE 5.5 DERIVED SCENE DEPENDENT 
ALGORITHM PARAMETERS FOR CZCS CHANNELS 1-4 

ALG. 1: Gordon's Epsilons (1...4) = (0.665, 0.886, 0.744, 1.00) 

ALG. 2: Smith and Wilson Epsilons (1...4) = (.0790, 0.972, 0.816, 1.00) 

ALG. 3: Aerosol Optical Depths, (1...4) = (0.327, 0.367, 0.350, 0.521) 

ALG. 4: Aerosol Vector Direction Cosines, V»(1...4) = 0.605, 0.540, 
0.433, 0.315) 

5.10 ALGORITHM COMPARISON AND PERFORMANCE ANALYSIS 

Three criteria were used to evaluate the performance of each 
atmospheric correction algorithm. First, derived or measured local 
optical depths were used with an atmospheric radiance model to estimate 
the aerosol portion of path radiance. These path radiances were 
compared with those predicted by each algorithm. Second, the upwelling 
CZCS spectra or that implied from known concentrations of constituents 
were compared with that obtained from each of the atmospheric correction 
algorithms. Third, the spatial patterns of haze known to be present in 
the raw uncorrected imagery were examined to see if they were removed by 
each algorithm. The residual patterns in the corrected imagery were 
also compared with known or suspected patterns of phytoplankton 
productivity and suspended particles. 
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Algorithm derived aerosol path radiances are compared with model 
estimates in Figure 5.6. The modified atmospheric water radiance model 
was used to estimate the aerosol path radiance over the Lake Superior 
calibration area. The atmospheric model values are shown here and they 
correspond closely to those produced by algorithm (4). The shape of the 
aerosol path radiance curve approximates theoretical considerations 
which suggest that aerosol path radiance is inversely proportional to 
wavelength. It is apparent from these results that both algorithms (1) 
and (2) severely underestimate the aerosol path radiance components. 
All four algorithms produced and equivalent estimate of the Rayleigh 
components of the satellite radiance. 

Upwelling radiance spectra were obtained for two test areas, one in 
southern Lake Michigan and the other in Lake Superior. These were 
compared with the spectra for each test site produced by the atmospheric 
correction algorithm. The upwelling radiance spectra used for the 
southern Lake Michigan test site were produced from water reflectance 
model using data collected coincidently with the CZCS's overflight. 

A comparison of the model predicted and algorithm predicted upwell- 
ing radiance for the southern Lake Michigan test site for algorithm (4) 
is presented in Table 5.6. Relative to Lake Superior these waters 

3 

contained high concentrations of chlorophyll-a pigments (4-29 mg/m ) and 
suspended minerals (0. 4-3.0 mg/1). Algorithm (3) produced similar 
results with the exception of band 4. Algorithm (1) produced the worst 
match to the model predictions. 
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TABLE 5.6. MODEL PREDICTED AND ALGORITHM PREDICTED 
UPWELLING RADIANCE (MW/CM 2 . SR . u) 


WAVELENGTH (NM) 


STATION* 

443 

520 

550 

670 

1 Model 

1.16 

1.39 

1.38 

0.40 

ALG. 4 

1.39 

1.58 

1.58 

0.49 

2 Model 

0.93 

1.07 

1.07 

0.25 

ALG. 4 

0.93 

0.93 

0.76 

0.19 

3 Model 

1.01 

1.16 

1.13 

0.25 

ALG. 4 

1.07 

0.92 

0.73 

0.23 

4 Model 

0.89 

0.81 

0.85 

0.25 

ALG. 4 

0.92 

0.88 

0.78 

0.20 

5 Model 

0.95 

1.22 

1.33 

0.57 

ALG. 4 

0.95 

1.31 

1.41 

0.43 


*A11 stations located near Grand Haven, Michigan 

Comparison of upwelling radiance predicted by each algorithm with 
that obtained from model calculations is shown in Figure 5.7 for the 
Lake Superior calibration area. Algorithms (2), (3), and (4) all 
approximate model values. Exceptions noted are that algorithm (2) 
produced a slightly higher value of Ly(x) at 443 nm and algorithm (4) a 
slightly higher value at 670 nm. The Gordon algorithm produced the 
poorest fit to the modeled radiances. 

The other means used for algorithm comparison and evaluation 
included transect and image analysis. A 50 km image transect was laid 
out in a west-east direction oriented perpendicular to the apparent haze 
pattern over Lake Superior. Figure 5.8 shows the uncorrected radiance 
at 443 nm along this transect with the corrected values generated by 
each algorithm. Algorithm (1), as applied, did not normalize the impact 
of aerosol haze. Each of the other three algorithms produced nearly the 
same result which appears to be free of the radiometric variations due 





to aerosol path radiance. Similar plots made for 520 nm and 550 nm pro- 
duced comparable results. Transect plots for 670 nm are shown in Figure 
5.9. The corrected radiance for algorithm (1) contains no upwelling 
contribution and clearly shows the effects of haze. The corresponding 
results for algorithm (2) and (3) imply haze normalization with only a 
very small upwelling radiance component. Algorithm (4) produced 
slightly higher radiance values implying a larger upwelling component. 
Transects selected for Lake Michigan waters produced similar comparisons 
between the four algorithms. Since these transect studies were made 
ovoer waters with minimum productivity these results suggest that the 
Gordon algorithm (Alg. 1) has little application for the Great Lakes 
region. 

Another type of spatial comparison of algorithm performance was 
made for algorithms (2), (3), and (4) by producing color-coded radiance 
maps in each band for both corrected and the original uncorrected data. 
A set of maps were made for both the Lake Superior and Lake Michigan 
test sites. These maps are shown respectively, as plates 1 and 2 of 
Appendix B. 

These maps illustrate the patterns of aerosols over each lake and 
the extent to which each algorithm was able to correct the subscene on 
an area wide basis. Examination of the uncorrected data indicated that 
large radiometric variations exist in both the aerosol and water com- 
ponents. The same water details were essentially absent in the 750 nm 
band where the variations due to aerosols were greatly attenuated. 
Algorithms (2) and (3) both used an aerosol radiance estimate in the 670 
nm band to estimate the corresponding radiance at the shorter wavelength 
bands. This process tended to eliminate spatial details in the upwell- 
ing radiance at 670 nm. An attempt was made to correct this situation 
by using the 750 nm band instead of 670 nm to estimate the aerosol 
component. However, the aerosol variations were represented by fewer 
signal counts at 750 nm than at 670 nm. The process caused stairstep 
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anomalies to appear in each of the first four CZCS bands. By compari- 
son, the approach used in the residual component algorithm (4) preserved 
much of the water radiance with Lake Superior. Algorithms (2), (3), and 
(4) all removed the haze patterns at 443, 520, and 550 nm. The residual 
patterns of upwelling radiance seem to be consistent with our expecta- 
tions, based upon suspected patterns of productivity and suspended 
sediments. As with the transect comparisons only algorithm (4) produced 
acceptable results at 670 nm. 

Results from these algorithms for Lake Superior displayed varying 
degrees of haze removal capability. Some residual haze features are 
apparent in the radiance corrected images produced by each algorithm. 
These features correspond to those portions of the subscene with the 
greatest concentration of haze. Algorithm corrected radiance compared 
well for portions of the east end of Lake Superior where the quantity of 
haze was the lowest. Over the entire lake residual variations due to 
haze were greatest for algorithm (3) and least for algorithm (4) which 
showed the best performance overall for the Lake Superior data set. 
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6.0 


DEVELOPMENT OF BIO-OPTICAL ALGORITHMS FOR THE GREAT LAKES 

A primary objective in this validation study is the development of 
working algorithms which are capable of transforming the satellite mea- 
sured radiance into surface concentrations of chlorophyll-a pigment and 
suspended sediment (minerals). In the previous section the development 
of atmospheric correction algorithms were discussed. This section 
focuses on the development of the companion water algorithms. These 
algorithms together are sufficient to convert spacecraft radiances in 
water constituent concentration values. The bio-optical algorithm is 
designed to transform the upwelling radiance values as estimated for a 
position just above the surface into concentration values for the 
selected biochemical constituents. The fundamental relationship 
developed in the bio-optical algorithm are based upon the absorption and 
scattering optical properties of each constituent. By contrast many 
previous algorithms were based upon empirical relationships between 
concentration values and remotely sensed radiances. This latter 
approach is usually valid for a limited range of environmental 
conditions since the optical properties are implicitly expressed in the 
algorithm. For example the NOAA chlorophyll algorithms derived for the 
open ocean are sensitive to small changes in chlorophyll concentrations 
in those open waters of the ocean, but have little application to 
coastal waters where a highly variable suspended particulates and 
dissolved organics can greatly influence and possibly dominate the 
observed radiance. 

The relationship between upwelling radiance and component concen- 
tration [c.] is usually expressed by the following equation. 

Log[c i ] = A : + A 2 Log(L u (x 1 )/L u (x 2 )) (28) 


69 


2p 

where and are determined by a least squares fit of a set of chlor- 
ophyll -a [C] or suspended sediment [SM] data and the radiometric data 
acquired at two appropriate wavelengths and X 2 » For CZCS a band 
1/band 3 ratio has been most commonly used where band 1 at 440 nm is 
highly sensitive to changes in pigment concentration and band 3 at 550 
nm is located near the hinge point of the chlorophyll absorption curves. 
In principle these algorithms cannot be applied without modification in 
waters having optical properties differing from the algorithm calibra- 
tion sites. The goal is to model the influence of these components on 
the reflectance in order to solve the inverse which would allow the 
derivation of component concentrations from reflectance. 

In order to systematically approach the development of new algori- 
thms it is necessary to understand the underlying radiative transfer 
processes and optical properties which influence the observed 
radiometry. The inherent optical properties of a passive medium include 
the absorption coefficient and the volume scattering function. From a 
practical standpoint it is necessary only to specify the total scatter- 
ing coefficient and the backscatter coefficient which is obtained by 
integrating the volume scattering function. For a water medium it is 
assumed that these inherent optical properties are a linear function of 
concentrations of various constituents. Thus, the total coefficient can 
be expressed as shown by the equation in section 5.4. 

There are several possible approaches to the determination of 
absorption and scattering cross sections for each of the constituents. 
First, one could isolate each constituent and measure the inherent 
optical properties for different concentrations and thereby estimate the 
cross section value. The principal problem with this approach is that 
the isolation may be difficult and could possibly alter the samples such 
that the measured conditions might not be representative of the natural 
waters. A second approach for determining optical cross sections is to 
measure the in situ apparent or inherent optical properties and 
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constituent concentrations for a large number of samples. Apparent 
properties can be transformed to inherent properties. Multiple linear 
regression can then be used to estimate the scattering cross sections 
for each component. However, biases could be introduced if there exists 
covariance between constituents. A third approach is to extract the 
optical cross sections indirectly from measured radiances and concentra- 
tion data. A radiative transfer model is used with numerical methods to 
interrelate the two sets of data and to extract the optical parameters. 

The advantage of this latter method is that it does not require 
extensive absorption and scattering measurements which are difficult to 
perform. The disadvantage is that the method may not converge to a 
correct or precise determination of the optical cross sections. Having 
prior knowledge of optical properties for one or more of the constitu- 
ents present can greatly improve chances of obtaining good approxima- 
tions for the others. Thus, while the specification of optical cross 
sections is recognzied to be a difficult task, one gains the necessary 
information to model the affect of constituent concentration on upwell- 
ing radiance and thereby, can gain a great deal of insight on water 
types throughout the Great Lakes. The empirical model, on the other 
hand, could not easily account for simultaneous changes in multiple 
parameters and may have little validity beyond the immediate area of 
surface truth collection. 

In principle the upwelling radiance at any point can be defined by 
the solution of the radiative transfer equation with proper boundary 
conditions and specified inherent optical properties. An analytical 
solution is not possible, but at least two approaches exist for 
obtaining an approximate solution. First, an exact numerical solution 
can be obtained using Monte Carlo and/or matrix operator methods [29]. 
Accuracy from these methods can be excellent but the costs make 
examination of a large number of solutions impractical. For the present 
study, the development of specific bio-optical algorithm for the Great 
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Lakes was based upon the Gordon power series model described previously 
in section 5.4. 

6.1 ANALYSIS OF OPTICAL PROPERTIES FOR THE GREAT LAKES 

The inherent optical data available to this project consisted of 
values reported in the literature, those produced from the Canadian Lake 
Ontario Study [16], and that measured by NASA/LaRC on the present 
program [10]. Previous to the present work and that for Lake Ontario 
several analyses [9, 20, 31] have been performed with marine waters to 
estimate the optical cross section for chlorophyll -a pigments using both 
laboratory and insitu optical measurements. Initially the greatest 
weight was logically given to the LaRC data since it was collected 
coincident with the 1980 sampling program. 

The optical measurements made by NASA/LaRC during the 1980 summer 
experiments were for lake samples. In this case the measured optical 
properties pertain to the particular mix of constituents in the lake 
sample. If sufficient number of measurements are made in this manner 
and if the principal constituents present are known then multiple 
regression can be used to derive the optical cross sections for a common 
parameter. While the present data are considered limited for this 
purpose a preliminary set of optical properties were derived for the 
Great Lakes. Of the twenty-one optical data set collected three Lake 
Erie samples contained sufficient quantities of sediment so as to 
saturate the optical measurement instrumentation. Samples collected 
from Green Bay in Lake Michigan, from western Lake Erie, and from 
western Lake Superior indicated that all these waters have somewhat 
distinctive local optical properties. Four of the six samples collected 
from Grand Haven area were essentially sediment free and the presence of 
low chlorophyll -a concentrations made absorption measurements difficult 
and results uncertain. For these reasons many of the original twenty- 
one samples were not considered suitable to be included in the 
regression analysis. 


2p 

In fact only eleven samples were used, nine from Lake Erie and two 
from Lake Michigan. Several linear regression models were formulated 
and tested against these selected optical measurements. These models 
utilized the available surface truth data including chlorophyll -a, 
phaeophytins, total residue, ashed residue, and volatiles [20]. Of 
these chlorophyll -a, chlorophyll-a plus phaeophytins (i.e., corrected 
values), total residue, and ashed residue were selected for regression 
analysis. Each regression model consists of four equations pertaining 
to the backscatter cross section and four equations pertaining to the 
absorption cross section. In each case the four equations correspond to 
the first four CZCS wavelengths (443, 520, 550, 670 nm). The optical 
model considered for these analyses describes the surface water mass to 
be a combination of unique organics as represented by chlorophyll-a and 
phaeophytin pigments and concentrations of unique inorganics as 
represented by the measurement of suspended minerals (SM). In addition 
each of the models contains a constant which includes absorption and 
s’cattering of pure water, dissolved organics, and possibly other 
unmeasured components. The two component model equations are written as 
follows. 


a = a 1 . [C] + a 2 . [SM] + a Q 
Bb = Bb x . [C] + Bb 2 . [SM] + Bb Q 


1 . 29 } 


Where [C] and [SM] are the concentrations of chlorophyll-a and suspended 
minerals, respectively and a and Bb are the absorption and backscatter 
cross sections. Statistics from these regressions are summarized in 
Table 6.1. Also included was a case where only the measurements for 
suspended minerals as ashed weight have been used. Corresponding 
optical cross sections derived from these regression analyses are given 
in Table 6.2. 
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The cross sections as derived for chlorophyll have the correct 
spectral shape but are only about one half the size of those derived for 
the Ontario Model and other reported results as summarized in Table 6.3. 
All of the samples included in this analysis were from very turbid 
waters with a high suspended sediment component which perhaps suffi- 
ciently dominated the measured optical properties so as to make these 
data less correlated to values of chlorophyll -a concentration. This 
dominance is confirmed by the single component regression results for 
suspended minerals (Table 6.2). Notice that the derived cross sections 
for suspended minerals are approximately the same as those obtained in 
the two component cases and where the multiple regression coefficient is 
0.999 and the standard error is greatly reduced. For this latter case 
the standard error of the estimate coefficient was approximately .02 of 
the coefficient value. Repetition of the two component model using 
total chlorophyl 1-a plus phaeophytin pigments produced almost the same 
cross sections for suspended minerals and slightly smaller values for 
chlorophyll-a than reported in Table 5.3. The large constant term in 
each regression equation also suggests that we have not properly 
accounted for the absorption due to chlorophyll-a and possibly other 
components as well. In any case the correlation with suspended sediment 
concentration is high and the cross sections are considered to be 
representative for the waters sampled. 


TABLE 6.3. CHLOROPHYLL ABSORPTION CROSS SECTIONS (m 2 /mg) 


Wavelength (nm) 

443 

520 

550 

670 

Smith and Baker [27] 

0.0293 

0.0143 

0.0090 

0.0113 

Morel and Prieur [9] 

0.0253 

0.0159 

0.0102 

0.0203 

Bukata, Jain et al [15] 

0.0354 

0.0240 

0.0173 

0.0086 

Present Work* 

0.0162 

0.0936 

0.0053 

0.0045 


*As derived from NASA/LaRC optical measurement. 
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The Canadian Lake Ontario optical water quality model was based 
upon a more comprehensive set of water quality parameters than available 
under the present program. This set allowed the specification of two 
additional components in the optical model for dissolved organics and 
non-living organics. The optical cross sections were not derived from 
measurements of the inherent properties as described above but rather 
were inferred from in situ measured apparent properties using the Gordon 
reflectance model [8]. The measured apparent properties included the 
diffuse reflectance, the irradiance attenuation, and the beam attenua- 
tion. Once the inherent properties of absorption and scattering were 
estimated they were deconvolved into their cross sections using 
regression techniques. 

In constructing an optical model for the present study the avail- 
able data were recognized to be insufficent to characterize any one of 
the Great Lakes. As discussed above the measured inherent optical data 
were insufficient to characterize local areas or single lakes, but alto- 
gether provided optical cross sections for suspended sediments. Further 
there was no complete set of optical data reported outside of that con- 
tained in the Lake Ontario model and it was not known how representative 
this model would be for all of the Great Lakes. Under these circum- 
stances a single set of optical cross sections was sought to represent 
the waters of the Great Lakes. It was conjectured that if the CZCS 
radiances can be predicted by a model at the positions where surface 
measurements of water parameters were collected, then the model optical 
cross sections are representative of those lake waters. The approach 
taken was to use a direct search method for unconstrained optimization 
where the water optical model which appeared to best fit the observed 
CZCS radiances could then be used for further investigation and 
algorithm development. 

Starting point for this modeling process was central Lake Superior 
where the chlorophyll -a concentrations are known to be very low. Typical 
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chlorophyl 1-a pigment levels for Lake Superior [27] average 1.1 mg/m 
while in the water radiance calibration area the level is typically 
lower at 0.8 mg/m . Chlorophyl 1-a concentrations are also very stable 
in this area during the spring to mid-summer period. Year to year pro- 
ductivity cycles in central Lake Superior are not expected to change 
dramatically based upon year to year measurements taken in similar 
waters of Lake Huron [20]. These latter samples represent the low range 
of chlorophyl 1-a and suspended mineral concentration values where as the 
measured boat values from the five CZCS dates represent mid to high con- 
centrations. Because of cloud cover open lake samples of CZCS radiance 
at satellite altitude could be extracted for Lakes Michigan, Superior, 
and Huron for the June 11 and July 23 images and also for Lake Superior 
for the August 6 scene for a total of six samples. The measured data 
available for this analysis consisted of the Lake Erie sites (July 14, 
2) and July 17, 2), Lake Michigan at Grand Haven (June 11, 5 and July 
23, 5), and Lake Superior (August 6, 5) for a total of 19 sampling 
stations. It was assumed that at each station there was a variable 
concentration of atmospheric aerosols and the aerosol type could also 
vary from scene date to scene date. The concentrations of chloro- 
phyl 1-a, suspended minerals, and dissolved organics were known or 
assumed. The unknown optical cross sections constitute eight unknowns 
and the unknown concentration of aerosol at each measurement site con- 
tributes a total of twenty five unknowns. Under these circumstances it 
is possible to write 25 x 4 non-linear differential equations in the 25 
+ 8 unknowns and any acceptable solution must be first bounded in the 
parameter space. One can start with an existing optical model or some 
notion of how the cross section spectra should appear although these 
cross sections would be expected to change if suspended sediments have 
different particle size distributions or absorption properties connected 
to differing parent soil materials (i.e., sand, silt, blue clay, and red 
clay). Further complexities to finding a solution are the uncertainties 
in the CZCS radiance measurements and possible insensitivities of 
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measured radiance to changes in constituent concentration. These latter 
insensitivities can occur at high concentrations of chlorophyll or 
suspended sediment. Although any attempt to obtain a general solution 
appears difficult an approximate solution is possible in the 
neighborhood of a suitable optical model. 

There are many well known approximate solution methodologies 
including LaGrange multipliers, Runge Kutta, and steepest descent. 
However, all of these methods require computation of derivative at each 
iteration. Since our solution to the radiative transfer equation is in 
the form of a computer model the computation of derivatives is a 
formidable task. An approach was sought which would find a solution 
without involving derivatives. The method selected was essentially that 
due to Hooke and Jeeves [11] which is a sequential technique where each 
step comprises two types of moves: EXPLORATORY and PATTERN. The first 

kind of move is designed to explore the local behavior of the objective 
function which can be written as 

F(a ij- Bb ij- c j> = G(a ij- B V c j> - L si 1 ■ i- 2 - 3 ’ 4 ; j * !’ 2 (3 °> 

where G is the model function and A., and B,. are the absorption and 

I vl * J 

backscattering coefficients for wavelength i and constituent j with 
concentration C- and L • is the satellite measured radiance. Introduc- 

J ** ' 

ing a starting point (a--, Bb^) we prescribe step lengths aA and ABb in 
each of the directions e^ , i =1 , 2, 3, 4. The exploratory stage is 
performed as follows. Set i=l and compute F(a^ + Aa, a 2 , a 3 , a^, Bb^ + 
ABb, Bb^, Bb 3 , Bb^) where Aa/ABb is constant. If F < F* accepts the 
trial point, set i=i+l, and repeat new trial point. Here F* is the 
value of the objective function from the previous base point (i.e., 
optical model). If F > F* reject trial point and compute F(a^ - Aa, a 2 , 
a 3 , a^; Bb 1 - ABb, Bb 2 , Bb 3 , Bb^). If again F < *F reject trial point, 
set i=i+l, and repeat from the beginning. After a series of exploratory 
moves we arrive at a new base point. If our objective function is less 
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at the base point then from the starting point F* then we substitute 
this new base point (i.e., pattern move) and restart the exploratory 
stage. Otherwise we return to the previous base point restart the 
exploratory moves reducing the size of Aa and ABb. These sequential 
moves are continued until satisfactory convergence is obtained. 

The solution methodology employed in the present study consisted of 
the following steps. 

1. Utilize principal component analysis to extract the spectral 
shape of the aerosol path radiance (i.e., as the first 
eigenvector) for each date and Great Lake. Use spectral shape 
factors to estimate atmosphere aerosol optical thicknesses 
produces the original extracted eigenvector. 

2. Set the initial optical cross section coefficients to those 
reported for the Lake Ontario five component model [16], Set 
absorption and backscatter coefficients for background waters 
to match oceanic clear water radiances as developed by Gordon 
[24]. 

3. Set model view geometry and solar geometry at each test pixel. 
Set concentration levels to expected or measured values. Iter- 
ate the atmospheric water radiance model over aerosol optical 
thickness until a best match is obtained between the observed 
and calculated CZCS band spectra at satellite altitude. Check 
the consistency and reasonableness of aerosol optical thicknes- 
ses derived for a sampling area (i.e., Lake Michigan, Grand 
Haven Stations). 

4. Investigate suitability of optical cross sections estimated 
from the NASA/LaRC optical measurements and as reported in the 
literature. 

5. Use the results from the above analyses to define the initial 
base point of the Great Lakes optical model. Apply the Hooke 
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and Jeeves exploratory and pattern search technique to converge 
on an acceptable optical model. 



6.2 DISCUSSION OF ANALYSES AND RESULTS 

Although this approach has a heuristic nature it is one which can 
produce a single Great Lakes optical model representative of the avail- 
able surface measurement and CZCS data sets. The optimality of this 
approach remains intuitive rather than that which can be substantiated 
easily by mathematical proof. While it is recognized that the derived 
model is based upon very little actual data for the size of the system 
and that a single model for the entire system is less appropriate than a 
series of models designed for each lake or sublake area the limitations 
of the present data set left little alternative. The extraction of 
needed optical parameters from the actual data sets by model search 
techniques is considered appropriate since there is little likelihood 
that extensive measurements of inherent optical properties can ever be 
obtained for a complex "coastal" system such as the Great Lakes. Thus, 
we view our approach to be consistent with the theoretical aspects of 
the water atmospheric radiative transfer process and to be practical in 
consideration of scarce optical and biological measurements. 

The inherent optical properties for the models considered are shown 
in Table 6.4 where model (I) is the Lake Ontario optical model, model 
(II) is the Lake Ontario model modified to include inherent absorption 
and scattering properties as measured in Lakes Erie and Michigan, and 
model (III) which is the result of five iteration steps from model (II) 
using the direct search methodology described above. Comparison of 
modelled and measured CZCS satellite radiances are made in Table 6.5 for 
each of the three models. A measure of the goodness of the fit of the 
measured and modeled data are the unbiased RMS errors computed for each 
of the CZCS wavelengths over the twenty-three stations and sites as 
shown in Table 6.6. The process of fitting the modeled to the measured 
data was done so initially to minimize the different at 443 nm. Because 
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TABLE 6.4 

OPTICAL CROSS SECTIONS FOR GREAT LAKES WATER RADIANCE MODELS 


MODEL I 


Wavelength(nm) a-C(m 2 /mg) a-SMfm^/mo) a-P0(in 2 /mq) a w 


1 wi. v 

443 

■ -,17 T'-'l/ 

0.0354 

— - L zai. 

0.0557 

“ 1 A 

0.0730 

0.020 

520 

0.0240 

0.0281 

0.0390 

0.028 

550 

0.0173 

0.0185 

0.0390 

0.037 

670 

0.0100 

0.0225 

0.0042 

0.370 


Bb-C(m 2 /mq) 

Bb-SM(.m 2 /mq) 

Bb-SM(m2/mq) 


443 

0.00199 

0.0328 

0.00 

0.0065 

520 

0.00182 

0.0474 

0.0 

0.00370 

550 

0.00241 

0.0525 

0.0 

0.00296 

670 

0.00175 

0.0333 

0.0 

0.0014 

MODEL II 

Wavelenqth(nm) 

a-C(m 2 /mq) 

a-SM(m 2 /mq) 

a-D0(m 2 /mq) 

!w 

443 

0.0354 

0.0764 

0.0730 

0.020 

520 

0.0240 

0.0636 

0.0390 

0.028 

550 

0.0173 

0.0577 

0.0390 

0.037 

670 

0.0100 

0.0556 

0.0042 

0.370 


Bb-C(m 2 /inq) 

Bb-SM(m 2 /mq) 

Bb-SM(m 2 /mq) 


443 

0.00199 

0.0312 

0.0 

0.0065 

520 

0.00182 

0.0284 

0.0 

0.00370 

550 

0.00241 

0.0287 

0.0 

0.00296 

670 

0.00241 

0.0250 

0.0 

0.0014 
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TABLE 6.4 (Cont'd.) 


MODEL III 
Wavelength 

2 

a-C(m /mq) 

a-SM(m 2 /mq) 

a-00(m 2 /mql 


443 

0.0350 

0.0764 

0.0730 

0.020 

520 

0.0240 

0.0636 

0.0390 

0.028 

550 

0.0170 

0.0577 

0.0390 

0.037 

670 

0.0280 

0.0556 

0.0042 

0.370 


Bb-C(m 2 /mq) 

Bb-SM(m 2 /mq) 

Bb-SM(m 2 /mq) 


443 

0.00200 

0.0315 

0.0 

0.0065 

520 

0.00190 

0.0285 

0.0 

0.00370 

550 

0.00180 

0.0240 

0.0 

0.00296 

670 

0.00175 

0.0200 

0.0 

0.0014 


Model I Lake Ontario Model 

Model II Lake Ontario Model modified to include suspended mineral 

absorption and backscatter. Cross sections measured for 
Lakes Erie & Michigan. 

Model III Modified Lake Ontario - Chlorophyll -a absorption and 

backscatter with measured absorption and backscatter for 
suspended mineral particles. 


C - Chlorophyll-a 
SM- Suspended Minerals 

a BB W - absorption and backscatter for background water 
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TABLE 6.5 COMPARISON OF MODELLED AND MEASURED 
CZCS SATELLITE RADIANCES 
2 

L(mW/sr era p) 


June 11 

. 1980 - 

Lake Mi chi 

qan Near Grand 

Haven 





A(nm) 

Model I 

Model II 

Model III 

CZCS 

Station 

1 

443 

6.73 

6.84 

6.84 

6.83 



520 

5.20 

4.60 

4.60 

4.71 



550 

5.04 

4.09 

3,84 

4.07 



670 

1 .81 

1 .73 

1 .57 

1 .64 

Station 

2 

443 

6.71 

6.67 

6.67 

6.67 



520 

4.70 

4.33 

4.33 

4.24 



550 

4.29 

3.77 

3.57 

3.45 



670 

1.57 

1 .51 

1 .43 

1 .45 

Station 

3 

443 

6.62 

6.56 

6.56 

6.61 



520 

4.73 

4.27 

4.27 

4.18 



550 

4.35 

3.71 

3.50 

3.36 



670 

1 .50 

1 .43 

1 .35 

1 .43 

Station 

4 

443 

6.68 

6.65 

6.65 

6.67 



520 

4.65 

4.31 

4.31 

4.24 



550 

4.26 

3.77 

3.56 

3.50 



670 

1.59 

1 .53 

1 .44 

1 .48 

Station 

5 

443 

6.71 

6.65 

6.85 

6.77 



520 

5.02 

4.45 

4.59 

4.71 



550 

4.95 

4.06 

3.91 

4.12 



670 

2.47 

2.16 

1 .82 

1 .70 

July 14, 

, 1980 - 

Lake Erie 

Near Sandusky 




Station 

1 

443 

6.63 

6.69 

6.69 

6.67 



520 

5.01 

4.60 

• 4.60 

4.76 



550 

4.76 

4.09 

3.86 

4.22 



670 

2.37 

2.19 

1 .89 

1 .61 

Station 

2 

443 

6.86 

6.86 

6.82 

6.83 



520 

5.17 

5.17 

4.69 

4.89 



550 

4.91 

4.91 

3.95 

4.41 



670 

2.53 

2.30 

1 .95 

1 .71 
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TABLE 6.5 ^Continued) 

August 6, 1980 (Continued) 



A(nm) 

Model I 

Model II 

Model III 

CZCS 

Station 4 

443 

6.05 

5.83 

5.83 

5.92 


520 

4.01 

3.67 

3.67 

3.83 


550 

3.28 

2.95 

2.82 

3.12 


670 

1 .33 

1 .21 

1 .16 

1 .15 

Station 5 

443 

6.14 

5.99 

5.98 

5.87 


520 

5.14 

4.03 

4.03 

3.67 


550 

4.75 

3.38 

3.18 

2.95 


670 

1 .36 

1.21 

1 .12 

1 .12 

July 23, 1980 - 

Central 

Lake Michiqan 




Station 1 

443 

7.36 

7.33 

7.33 

7.33 


520 

4.67 

4.44 

4.44 

4.40 


550 

3.79 

3.52 

3.43 

3.37 


670 

1.24 

1.23 

1.22 

1.26 

July 23, 1980 - 

Central 

Lake Huron 




Station 1 

443 

6.74 

6.70 

6.70 

6.62 


520 

4.29 

4.04 

4.04 

4.02 


550 

3.49 

3.22 

3.12 

3.16 


670 

1 .09 

1.07 

1 .06 

1 .32 

July 23, 1980 - 

Central 

Lake Superior 




Station 1 

443 

7.03 

7.09 

7.09 

7.04 


520 

4.44 

4,29 

4.29 

4.25 


550 

3.60 

3.41 

3.31 

3.24 


670 

1 .16 

1 .19 

1 .17 

1 .32 

June 11. 1980 - 

Central 

Lake Michiqan 




Station 1 

443 

6.82 

6.79 

6.79 

6.85 


520 

4.43 

4.22 

4.21 

4.31 


550 

3.70 

3.44 

3.32 

3.46 


670 

1 .30 

1 .29 

1 .27 

1 .49 

June 11, 1980 - 

Central 

Lake Huron 




Station 1 

443 

6.43 

6.40 

6.40 

6.56 


520 

4.38 

4.15 

4.14 

4.11 


550 

3.64 

3.37 

3.27 

3.26 


670 

1 .42 

1 .40 

1.39 

1 .28 

June 11 , 1 980 - 

Central 

Lake Superior 




Station 1 

443 

7.18 

7.15 

7.15 

7.14 


520 

4.96 

4.72 

4.72 

4.55 


550 

4.20 

3.93 

3.83 

3.64 


670 

1 .73 

1 .71 

1 .70 

1 .78 
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TABLE 6.5 (Continued) 


July 17, 1980 - Lake Erie Near Sandusky 



A(nm) 

Model I 

Model II 

Model III 

CZCS 

Station 1 

443 

8.48 

8.36 

8.36 

8.38 


520 

7.41 

5.81 

5.81 

6.77 


550 

7.56 

5.18 

4.87 

5.67 


670 

4.43 

3.20 

2.58 

2.69 

Station 2 

443 

7.82 

7.89 

7.88 

7.84 


520 

5.60 

5.16 

5.16 

5.72 


550 

5.24 

4.52 

4.29 

4.71 


670 

2.66 

2.43 

2,08 

1 .80 


July 23, 1980 - Lake Michigan Near Grand Haven 


Station 1 

443 

7.15 

7.17 

7.12 

7.15 


520 

4.78 

4.47 

4.47 

4.38 


550 

4.32 

3.87 

3.67 

3.57 


670 

1.59 

1 .52 

1 .42 

1 .37 

Station 2 

443 

7.13 

7.11 

7.10 

7.10 


520 

4.67 

4.42 

4.42 

4.34 


550 

4.12 

3.78 

3.60 

3.52 


670 

1 .49 

1 .44 

1 .37 

1 .32 

Station 3 

443 

7.18 

7.18 

7.17 

7.15 


520 

4.47 

4.44 

4.40 

4.41 


550 

3.97 

3,92 

3.70 

3.64 


670 

2.10 

2.08 

1 .68 

1 .41 

Auqust 6, 1980 - 

Lake 

Superior Near 

Duluth 



Station I 

443 

6.03 

5.96 

5.96 

5.92 


520 

4.64 

3.96 

3.96 

3.92 


550 

4.50 

3.49 

3.25 

3.28 


670 

1 .89 

1.64 

1 .38 

1 .31 

Station 2 

443 

6.10 

5.96 

5.96 

5.87 


520 

5.06 

4.00 

4.00 

3.73 


550 

4.69 

3.37 

3.17 

3.02 


670 

1 .36 

1 .21 

1 .12 

1 .09 

Station 3 

443 

6.02 

5.98 

5.98 

5.97 


520 

4.26 

3.88 

3,88 

3.99 


550 

3.80 

3.26 

3.07 

3.31 


670 

1 .49 

1 .40 

1 .29 

1 .38 
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TABLE 6.6. RMS ERRORS FOR MODEL PREDICTED RADIANCE* 


Wavelength (nm) 

Model I 

Model II 

Model III 

443 

0.0988 

0.0574 

0.596 

520 

0.413 

0.272 

0.267 

550 

0.851 

0.265 

0.265 

670 

0.497 

0.319 

0.156 


♦Calculated RMS errors based upon results reported in Table 6.5 


of the larger radiance values at this wavelength and the sensitivity of 
any subsequent predictive technique to the radiances in this band. The 
initial adjustments were made by varying the atmospheric aerosol optical 
thickness at each measurement site. An estimate of the aerosol optical 
thickness for each scene date can be obtained by subtracting the known 
Rayleigh and ozone optical thicknesses from the measured total value as 
reported in Table 3.6 of Section 3.6. These estimates are compared in 
Table 6.7 with model derived values of the aerosol path radiance as 
obtained with the final optical model. As would be expected the model 
predicted station to station variations in the aerosol optical thickness 
were very small or nil for each of the primary sampling sites. The RMS 
errors associated with the Lake Ontario model I are large especially at 
CZCS wavelengths 520, 550, and 670 nm. The inclusion of optical cross 
sections for suspended sediments derived from the NASA/LaRC data as 
model (II) represented a substantial reduction in the RMS errors and 
especially in bands 1 and 3. The least reduction was made in the 670 nm 
band. Between models (II) and (III), four exploratory moves were made 
at each wavelength in attempt to reduce the RMS errors. Little 
reduction was obtained except for the 670 nm band where RMS errors were 
reduced further by about 50% in model (III). Attempts to further reduce 
these errors in band 2 and 3 were fruitless and suggest that major 
differences in optical properties exist at these wavelengths throughout 
the Great Lakes. These differences seemed most apparent with July 14 
and 17 Lake Erie sampling sites, the radiance values for which implied 



TABLE 6.7 

COMPARISON OF MODEL PREDICTED AND MEASUREMENT DERIVED AEROSOL PATH RADIANCE 

? 

L w (nw/cm> -sr-p) 


Wavelength(nm) 


Date 

Site Location 

443 

520 

550 

620 

June 11 

Grand Haven, Lake Michigan 
Model Predicted 
Measurement Derived 

1.76 

1.33 

1.13 

0.840 



0.723 

0.557 

0.465 

0.341 

July 14 

Catawba, Lake Erie 
Model Predicted 
Measurement Derived 

2.06 

2.38 

1.54 

1.26 

1.26 

1.44 

0.832 

0.958 

July 17 

Catawba, Lake Erie 
Model Predicted 
Measurement Derived 

1.29 

1.89 

0.895 

1.30 

0.754 

1.10 

0.567 

0.831 

July 23 

Grand Haven, Lake Michigan 
Model Predicted 
Measurement Derived 

0.366 

1.94 

0.274 

1.44 

0.230 

1.22 

0.157 

0.830 

Aug 6 

Duluth, Lake Superior 
Model Predicted 
Measurement Derived 

0.212 

1.40 

0.191 

1.26 

0.106 

1.09 

0.038 

0.833 
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an exceptionally large backscatter component at 550 nm. Adjusting to 
these samples would throw all the others off greatly increasing the RMS 
errors for these bands. Thus, optical properties of model (III) were 
considered for all practical purposes to have the greatest "overall 
representativeness" of the optical properties of the Lakes with little 
opportunity to improve the result using this direct solution search 
methodology. 

6.3 USE OF THE GREAT LAKES OPTICAL MODEL TO DERIVE A BIO-OPTICAL 

ALGORITHM 

The optical cross sections derived above provided a means to simu- 
late the upwelling radiance from Great Lakes waters with a wide variety 
of component concentrations. Further the irradiance reflectance can be 
readily calculated for each CZCS band as a function of the concentration 
of chlorophyl 1-a pigments and suspended sediments. The spectral char- 
acteristics of the irradiance reflectance function can be depicted with 
iso-chloropleths for each pair of wavelengths. Figure 6.1 and 6.2 show 
calculated subsurface reflectance for the Great Lakes Optical model 
(III). Dissolved organics are assumed, for this case, to be 2.0 rng/t. 
Each of the panel figures has nine curves of increasing suspended 
mineral concentration and constant value of chlorophyll pigment concen- 
tration (0.0, 1.0, 2,0, 5.0, 10.0, 20.0, and 1000.0 yg/t). Each figure 
also contains four curves of increasing chlorophyll with constant values 
of sediment (0.0, 2.0, 10.0, and 50.0 mg/£). Each panel figure is a 
projection of the reflectance space and together suggest the reflectance 
space is a three dimensional hyperplane which is nearly perpendicular to 
the 520/550 plane as shown in Figure 6.2 as panel 3. This result 
suggests that there is little bio-optical information in this latter 
band pair which is contrary to the results obtained for the clear ocean 
where these bands were successfull used to predict chi orophyl 1 -a 
concentration [20]. 
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«!550n.] (I) »(4«Jn.) (I) 



SUBSURFACE IRRADIANCE REFLECTANCE MODEL III 



*(670n») (I) (») 

(c) (d) 


Figure 6.1 Subsurface irradiance reflectance (percent) at CZCS wavelengths 
(443 nm, 520 nm,' 550 nm , 670 nm) as predicted by the Great Lakes 
Optical model. Each figure has eight parametric curves of 
increasing suspended mineral concentration (0.0-50.0mn/n) with 
constant values of chlorophyll pigment concentration (0.0, 1.0, 

2.0, 5.0, 10.0, 20.0, or 1000 yg /i). Each figure also contains 
four parametric curves of increasing chlorophyll a (0.0-1000. Og/t) 
at constant values of suspended mineral concentration (0.0, 2.0, 

10.0, 50.0 mg/t). 
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SUBSURFACE IRRADIANCE REFLECTANCE MODEL III (Continued) 



K(WOi») (1) 


(b) 

Figure 6.2 Subsurface irradiance reflectance (percent) at CZCS wavelengths 
(443 nm, 520 nm, 550 nm, 670 nm) as predicted by the Great Lakes 
Optical model. Each figure has eight parametric curves of 
increasing suspended mineral concentration (0.0-50.0 mg/ it) with 
constant values of chlorophyll pigment concentration (0.0, 1.0, 
2.0, 5.0, 10.0, 20.0, 50.0, 1000 ug/t) at constant values of 
suspended mineral concentration (0.0, 2.0, 10.0, 50.0 mg/t). 
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The orientation of the band pair projection and the corresponding 
spacing between iso-chl oropleths curves suggests the potential of that 
individual band pair to resolve various constituent concentration 
levels. For example the (1,3) band pair (panel 2, Figure 6.1) provides 
the greatest resolution of chlorophyll and suspended sediment 
concentrations in the concentration range of 0 to 10 pg /l chlorophyll 
and 1 to 10 mg/t of suspended sediments. For greater concentrations of 
suspended sediments it appears that resolution of the chlorophyll con- 
centration would become more difficult especially for low concentration 
values. It is apparent from the spatial features of the subsurface 
irradiance reflectance function that each band pair coordinate repre- 
sents a unique concentration of chlorophyll -a and suspended sediments. 
If the concentration of dissolved organics were also involved then a 
triple in 3-space would be necessary to describe a set of unique 
concentrations. While the reflectance model provides a one to one 
mapping from component concentration space to band pair reflectance 
space the inverse mapping function is needed for purposes of extracting 
component concentrations. The complexity of the reflectance model makes 
any analytical form of the inverse function an impractical for purposes 
of extracting constituent concentrations. Since the model is 
essentially a numerical approximation to Monte Carlo calculations a 
numerical approximation to the inverse function seems appropriate. 

Points in the two band reflectance space (1,3) representing consti- 
tuent concentrations can be located with a vector of length equal to 
( ) + (R 2 ) and orientation equal to tan (R^/R^. Thus, two spectral 
variables, i.e., two band brightness and band ratio, are needed to 
resolve the two concentration components. The strong non-linearities in 
mapping reflectance to concentration values suggests that approximation 
functions are needed over subregions of the reflectance space. The 
approach taken was to partition the space using the brightness variable 
and then obtain polynomial approximations in the band ratio variable to 
the desired concentration pair. The approximation polynomials were 
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derived as follows first the Great Lakes optical model III was used to 
simulate the irradiance reflectance such that the simulation data set 
(3000 samples) was representative of the two band projection space and 
constituent combinations found in Great Lakes waters. These data were 
then categorized into ten separate regions as level ranges of the band 
pair reflectance brightness (i.e., two band vector modulus). 
Multivariate regression techniques were then used to estimate 
constituent concentration polynomial coefficients using the (1,3) band 
ratio as the independent variable. Polynomial coefficients were derived 
for both suspended sediments and chlorophyll-a pigments. A third 
degree, three term polynomial was found to be sufficient to describe 
constituent concentrations. Figure 6.3 shows an example of radiance 
versus concentration curves for chlorophyll-a pigment. Each of the 
concentration polynomials has the following form. 

C(R 1 /R 2 ) = a R ] (l,3) + a 2 R 2 (l,3) + a 3 R 3 (l,3) ( 31 ) 

Where R ( 1 ,3 ) is the band 1/band 3 irradiance reflectance ratio. Each of 
the curves essentially applies to a separate range of suspended sediment 
concentration. Note that for small concentrations of chlorophyll the 
curves of Figure 6.3 approximately overlay one another since the concen- 
tration of suspended sediment can be expected to be very small. At the 
other extreme the curves are highly separated based on the influence 
suspended sediments. As these polynomials curves become more vertical 
the ability of the band ratio to resolve chlorophyll concentration 
becomes greatly diminished as was apparent with the plots shown in 
Figure 6.1. 

Polynomials were also derived for a band 1/band 4 ratio algorithm 
which also appears from the simulation analysis to have capabilities of 
identifying chlorophyll suspended sediment constituent concentration 
pairs. Examination of simulation data sets for band ratios (1,2), 
(2,4), and (3,4) all were shown to have much less capability to extract 
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the desired concentrations levels, as indicated from the band pair 
reflectance space projections. The polynomials generated above were 
incorporated into a b.io-optical algorithm as described below. 


6.4 CZCS IMAGE PROCESSING FOR CHLOROPHYLL-A AND SUSPENDED SEDIMENT 

CONCENTRATIONS 

The five CZCS scenes were processed for two types of output. First 
a 3 x 3 array of pixels for each sampling site was processed on a pixel 
by pixel basis and the results averaged to produce a CZCS extracted 
concentration value for chlorophyll -a pigments and suspended sediments. 
Comparisons were then made to the available surface truth data. Second, 
a series of twelve subscenes were processed with individual scene, solar 
geometrical, and atmospheric aerosol characteristic parameters (as made 
available from pass 1 or by surface measurements). The approximate 

coverage of each subscene is shown in Figure 6.4 and subscene processing 
parameters are listed in Table 6.8. Processed subscenes were subsequen- 
tly used to create color coded level slice constituent concentration 

•> 

maps for chi orophyl 1 -a pigment, < 1.0 to 40 mg/m , and suspended sedi- 
ments, < 0.1 to 20 rng/t. The overall processing scheme is outlined with 
the data flow diagram shown as Figure 6.5. Pixel by pixel processing 
was performed using the ERIM QLINE image processing software with major 
sequential processing steps handled by individual image processing 
modules. Modules developed specifically for this CZCS validation study 
included: 

1. SPLIT - Cloud/land/water discrimination module. Option to 
black out land and cloud features. 

2. ATMOS - Atmospheric correction routine which removes the 
spatially varying aerosol path radiance and view angle effects. 
Pass 1 calculates proper byte to radiance conversion factors 
and aerosol characteristic parameters. Pass 2 is used to 
calculate upwelling water radiance or nadir radiance normalized 
and corrected to the top of a Rayleigh atmosphere. 
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FIGURE 6.4 CZCS GREAT LAKES SUBSCENES 
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CZCS LEVEL I DATA TAPES 


c 


QLINE DATA REFORMATTING AND SUBSCENE 
SELECTION 


) 


ground point control 
mapping polynomial generation 

NEAREST NEIGHBOR IMAGE RESAMPLING 



FIGURE 6.5 


OIAGRAM OF IMAGE PROCESSING SCENE 








3. WATPRM - Uses output from ATMOS to estimate concentration of 
chlorophyll -a pigments and suspended sediments. Options for 
bands (1,4) and (1,3) estimates and selected outputs as four 
channels. 

Input CCTs were geometrically corrected and resampled at 500 meters 
in a polyconic projection. Output tapes contained image files of 
chlorophyl 1-a and suspended sediment concentration values. Data flow 
diagrams for each of the above modules are shown as Figures 6.6, 6.7, 
and 6.8. 

6.5 ANALYSIS AND INTERPRETATION OF RESULTS 

Measured and algorithm predicted concentrations for chlorophyll-a 
pigments and suspended sediments for each of the sampling stations are 
listed in Table 6.9. Results are shown for both the (1,3) and (1,4) 
band ratios. While in many cases the predicted concentrations., from 
these ratios are similar the (1,3) ratio was felt to give better results 
and have a wider application to Great Lakes waters and especially in 
clear waters where the radiance in 670 nm band is very low and 
unreliable for extraction of concentration values. The average accuracy 
of NOAA type ocean CZCS chlorophyll algorithms have been reported to be 
within 0.5*1 og(C) [7]. The average log(C) coefficients for the present 
study were determined to be: 


June 11, 1980 
July 14, 1980 
July 17, 1980 
August 6, 1980 


Grand Haven, Lake Michigan 
Western Lake Erie 
Western Lake Erie 
Duluth, Lake Superior 


0.51LOG(C) 

0.13LOG(C) 

0.68LOG(C) 

0.53L0G(C) 


While the above results are based upon a small number of stations 
where surface truth data were gathered, they suggest that these Great 
Lakes CZCS algorithms may have comparable performance characteristics to 
that reported for the ocean experience. 

The results shown in Table 6.9 for July 23, 1980, Grand Haven, Lake 
Michigan are considered to be anomalous. Examination of the CZCS data 



Function: Separate cloud haze and land pixels 

Required Inputs: Faw count values for Bands 5 and 6 



FIGURE 6.6 DIAGRAM OF SPLIT QL1NE MODULE 
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Function: 


Pass 1: Extraction of Scene Dependent Atmospheric Parameters 

Pass 2: Remove Spatial Haze and View Dependent Components 



FIGURE 6.7 DIAGRAM OF ATMOS : ALGORITHM 4 QLINE MODULE 
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Function: Convert ATMOS corrected data to estimate of Chlorophyll -a and Suspended Mineral 


Concentrations 



FIGURE 6.8 DIAGRAM OF WATPRM QlINt MODULE 
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TABLE 6.9 



COMPARISON OF MEASURED AND ALGORITHM PREDICTED 
CONCENTRATIONS FOR CHLOROPHYLL PIGMENTS 
ANO SUSPENDED SEDIMENTS 
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P.31 R41 


R31 

R41 

June 11 

1 

6.91 

22.0 7.40 

1 .09 

05 

T75 

(Lake 







Michigan) 

2 

5.24 

3.2 3.2 

0.40 

0.8 

0.6 


3 

4.55 

3.0 3.8 

0.45 

0.8 

0.8 


4 

5.99 

3.4 4.4 

0.41 

<0.2 

0.8 


5 

29.3 

30.0 11.0 

2.95 

5.0 

1 .4 

July 14 

1 

18.1 

18.8 10.8 

2.01 

4.6 

2.2 

(Lake Erie) 








2 

21.1 

22.6 11.2 

2.2 

6.8 

3.0 

July 17 

1 

43.0 

Sat & 40 

21 .7 

Sat & 20 


(Lake Erie) 








2 

22.9 

14.6 5.2 

2.73 

2.6 

1 .0 

July 23 

4 

7.7 

1.8 2.0 

0.52 

0.8 

0.4 

(Lake 







Michigan) 

5 

6.2 

1.8 1.8 

0.31 

0.8 

0.4 


6 

78.5 

2.2 2.0 

0.42 

0.8 

0.4 

August 6 

1 

13.1 

6.6 11.2 

2.05 

0.2 

0.6 

(Lake 







Superior) 

2 

1 .72 

3.8 2.6 

1 .0 

<0.2 

0.8 


3 

7.16 

6.8 11.8 

0.67 

0.2 

0.8 


4 

3.88 

3.8 2.4 

0.2 

0.2 

<0.2 


5 

1 .43 

3.2 4.2 

1 .0 

<0.2 

<0.2 

R31 = Lu(550/Lu(443) 

R41 

= Lu(670)/Lu(443) 
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in the immediate area of each sampling site showed both the count values 
and band to band spectral shape to be consistent with low concentrations 
of chlorophyll and suspended sediments. The measured concentrations on 
the other hand were medium to very high. No sizable errors in the image 
registration and geometry could be established. Spectral shapes of the 
subsurface radiometer measurements are consistent with the CZCS derived 
shaped for the upwelling radiance and do not indicate the presence of 
large chlorophyll or sediment concentrations. Since all of the Grand 
Haven sampling stations were located in the Grand River plume or in 
coastal waters influenced by the river it is possible that the optical 
properties of the river waters substantially between the June 11 and 
July 23 sampling dates. 

Water samples collected at station 1, July 17, 1980, Western Lake 
Erie showed very high concentrations of chlorophyll -a and suspended 
sediments. CZCS derived radiances were either saturated or implied 
spectral brightness values were larger than the upper bounds set in the 
WATPRM module. As a result no reliable estimate could be extracted from 
the CZCS data for station 1. 

Color coded images for each processed subscene depicting either 
chlorophyl 1-a pigment or suspended sediment concentrations are shown in 
Appendix B. A total of twelve different subscenes were processed using 
the subscene parameters listed in Table 6.8. A brief description and 
interpretation of selected concentration maps on a lake by lake basis 
follows. 

Concentration maps are grouped as whole lake and sublake scenes. A 
total of fifteen colors were used to code a very wide range of CZCS 
extracted chlorophyll and suspended sediments concentrations. The color 
codes are graduated to show very small changes at the low end and very 
large changes at the high end. This scheme was considered consistent 
with the 0.5 Log(C) errors and brings out both open lake and nearshore 
features in the concentration of chlorophyl 1-a pigments and suspended 
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sediments. These color codes were applied in exactly the same way to 
the output of the WATPRM module for each of the scene dates. 

Lake Superior. August 6. 1980 . The Keweenaw current is clearly 
visible in the whole lake subscene. Also there are subtle differences 
in the chlorophyll concentration ( . 1- . 2 pg/t) and primary productivity 
from the western end of the lake as opposed to the central region just 
north and east of the Keweenaw Peninsula. These waters may be at the 
center of a large clockwise circulation pattern. The thermal maps 
consistently show these central lake waters to be cooler. It has been 
conjectured by some that these waters never stratify throughout the ice 
free period and, therefore, probably do not mix extensively with the 
more productive and nutrient rich coastal waters. The western Lake 
Superior subscene shows a large concentration gradient in the vicinity 
of the Duluth harbor and river discharge. Also visible are less 
productive colder waters upwelling along the Minnesota shoreline north 
of Duluth. Black Bay on the Canadian North Shore stands out as highly 
productive. Since this bay is very shallow, some of the radiance 
associated with productivity may in fact, be bottom reflected energy. 

Lake Michigan, June 11, 1980 . The structure of the nearshore 
thermal bar system is quite evident in these maps with high concentra- 
tions of chlorophyll and suspended sediments essentially ringing the 
lower basin shoreline. Since each map was generated with processing 
using a single set of atmospheric and solar geometrical parameters small 
differences can be observed in the resulting products processed from the 
same CZCS scene. Case in point are the whole lake and southern Lake 
Michigan concentration maps, the latter of which shows slightly higher 
concentrations than the whole lake map for the central lake waters. 

Chlorophyl 1-a pigment concentration values representative of the 
open lake waters, bays, and other sublake regions are summarized in 
Table 6.10 along with corresponding values reported in the 1976 Upper 
Great Lakes Reference Group Study [12]. 
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TABLE 6.10 

COMPARISON OF CZCS DERIVED CHLOROPHYLL-A (mg/m 3 ) 
PIGMENT CONCENTRATIONS WITH IJC SUMMARY RESULTS 



Predicted 

by CZCS 

I JC/ULRG (12) 

Location 

Summer, 

1980 

Circa. 1976 

Duluth Area 

2.0 - 

3.0 

2.0 

Lake Superior 

1.1 - 

1.2 

0.95 

N. Lake Huron 

1.2 - 

1.4 

1.4 

North Channel 

1.2 - 

1.4 

1.7 

Georgian Bay 

1.1 


1.2 

Saginaw Bay 

15.0 - 

30.0 

15.7 

S. Lake Huron 

1.4 


1.4 

W. Lake Erie 

15.0 - 

20.0 

11.1 

C. Lake Erie 

1.7 - 

3.0 

5.5 

E. Lake Erie 

1.7 - 

3.0 

4.3 

Lake Ontario 

5.0 


4.8 

Lake Michigan 

1.4 - 

1.7 

1.3 


IJC ULRG International Joint Commission Upper Lakes 
Reference Group 
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These latter values represent summaries of ship gathered chloro- 
phyll data over a several year period of the early seventies. The 
reported values which can be considered typical of those water bodies 
are in general remarkably close to the values derived from the CZCS data 
and especially in the clear low productive waters of the Great Lakes 
where the summer productivity levels would be expected to change very 
slowly if at all. In highly productive areas such as Saginaw Bay and 
the western basin of Lake Erie the reporting of a typical value is more 
dubious as is the comparison with previous data. The productivity of 
these latter waters can be expected to change dramatically over a short 
time. 
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7.0 

CONCLUSIONS 

This program to validate the use of CZCS for application to Great 
Lakes waters has resulted in the development of quantitative algorithms 
designed to correct the radiometric and geometric distortions which are 
inherent in the CZCS data sets as made available from NOAA. This study 
which has been conducted in two separate phases as a data collection and 
data analysis and algorithm development phase has led to the following 
set of conclusions and recommendations. 

7.1 GENERAL CONCLUSION 

The prospects for using satellite data and specifically that from 
the CZCS to study the Great Lakes ecosystem have been greatly enhanced 
by the advancements in quantitative methods developed under this study. 

Changing prospects to realities requires that satellite data be brought 
to bear on some major problem in the region. The Great Lakes represent 
large and diverse ecosystems with enormous spatial heterogeneity. 

Because of the heterogeneity in water and atmospheric optical properties 
over the Great Lakes, extraction of meaningful chlorophyll and suspended 
sediment information is considerably more difficult there than in the 
deep ocean. It is the conclusion of this study that algorithms can be 
developed to adequately handle the conditions present in the Great Lakes 
environment although these algorithms will be generally more complex 
than those which have been developed for the ocean case. Because of 
limited experimental data the approach taken was to develop a single 
algorithm for the entire Great Lakes waters. Predicted chlorophyll-a 
pigment concentrations in the open lake waters of Huron, Michigan, and 
Superior were comparable to long term average values. On the other 
hand, concentrations of chlorophyll and suspended mineral were more 
difficult to predict for the highly turbid and productive waters of Lake 
Erie, Saginaw Bay, Western Lake Superior. Much of the surface truth 
data collected on this program was collected in these highly turbid 
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areas which may have led to large CZCS sampling errors. Most CZCS 
predicted chlorophyll pigment concentrations were within 0.5 Log n [C] of 
the measured values which is consistent with the experience of other 
investigators [5,6,7], There are not sufficient data to establish a 
separate optical model for each lake or sublake region but in situ 
bio-optical measurements suggested that such a modelling approach would 
be appropriate. Largest errors in the predicted concentration were 
found at sampling sites located in or near sediment plumes from 
tributary discharge. These errors may have resulted from simultaneous 
sampling of multiple water masses. In this case the CZCS instantaneous 
field of view is ultimately averaging the upwelling radiance from the 

plume with that from adjacent waters having less turbidity. Thus 

inaccuracies in geographical location of any pixel and the large IFOV 
imply a need for extensive surface sampling where water quality 

gradients are significant. 

7.2 SPECIFIC CONCLUSIONS 

1. A geometric correction algorithm was developed based upon a 
CZCS scanner model and was found to be capable of resampling 
CZCS pixels to within approximately plus or minus two pixels 
RMS. 

2. A comparison of existing atmospheric correction algorithms with 

one developed for the Great Lakes showed (a) that significant 
quantities of radiance are upwelled in the 670 nm band invali- 
dating the possible use of the NOAA algorithms developed for 
ocean applications, (b) that the residual component method, 

developed for this project, was able to remove the radiometric 
effects of haze when those effects are severe. However, under 
those same severe conditions errors in estimating the upwelling 
radiance will be larger. 

Under conditions of light haze. Smith and Wilson [5] iterative 
algorithms, the pseudo optical depth algorithm [6], and the 
residual component algorithm all appeared capable of removing 
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haze in CZCS bands 1, 2, and 3. The latter algorithm appears 
to be able to normalize the severe haze conditions in band 4 as 
wel 1 . 


3. A principal component technique was used to extract the neces- 
sary scene dependent aerosol optical parameters for each algor- 
ithm from the CZCS data minimizing the need for independent 
ground measurements. The ability of this technique to work is 
dependent on two basic requirements. First, there must be some 
water masses in the scene which have low and stable levels of 
chlorophyll-a pigment concentration. Further, it must be 
possible to estimate the upwelling radiance for the calibration 
areas under a variety of illumination conditions. This need 
can best be satisfied with a combination of upwelling radiance 
measurements and/or reflectance model calculations. Second, 
the scene must contain variable concentrations of aerosol haze 
over the required calibration areas in order to extract the 
aerosol characteristic vector. 

4. Chlorophyll and sediment prediction algorithms were based on a 
reflectance model developed for the Great Lakes. A two para- 
meter algorithm was found to be necessary to specify both the 
chlorophyll and suspended mineral concentration from the same 
extracted upwelling radiance values. These parameters are a 
CZCS band one to band three ratio and the normalized spectral 
brightness in bands one and three. 

5. Chlorophyll absorption coefficients derived from the NASA/LaRC 
optical data set were found to be approximately one half of 
those obtained from the Lake Ontario study. Further, the 
backscatter coefficients were only about one tenth those 
obtained by the same study. Chlorophyll absorption and 
scattering cross sections as obtained by this regression 
analysis were considerably less than those reported in the 
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literature [9, 16]. However, most samples included in this 
analysis were collected from turbid waters where the optical 
properties are dominated by suspended solids. 

6. Simulation studies which described the reflectance space for 
the Great Lakes optical model suggest that concentrations of 
chlorophyll -a and suspended sediment are highly correlated to 
ratios of bands one to four or one to three. A band two to 
three ratio, on the other hand, offers little sensitivity to 
the concentration of chlorophyll. This result is somewhat 
contrary to the marine experiernce where Gordon has suggested 
that this ratio will prepdict such concentrations in the 1.5 to 
20 pg/fi. range. 

Lakes Ontario and Erie as well as Saginaw Bay and Green Bay have 
high concentrations of dissolved organics (approximate 2.5 mg/fc). In 
these water bodies the influence of dissolved organics on volume 
reflectance is minimal except for extreme cases. In many parts of Lakes 
Huron, Superior and Michigan the influence of dissolved organics is 
pronounced since the concentrations of chlorophyll and suspended 
sediments are low. For example in central Lake Superior changes in the 
dissolved organic concentration on the order of 0.2 mg/fc would have a 
very noticeable affect on the upwelling radiance spectra. This 
conclusion is based upon reflectance model studies and interpretations 
made of CZCS spectral signatures found in the open clear waters of the 
upper Great Lakes. 

7.3 RECOMMENDATIONS 

This study and our assessment of the requirements for Great Lakes 
monitoring has led to making the following recommendations. 

1. In order to refine CZCS algorithms for the Great Lakes a 
greater body of data are needed on the optical properties of 
Great Lakes waters. If possible in situ measurements should be 
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obtained of necessary apparent and inherent absorption and 
scattering properties. Accuracies of CZCS extracted 
concentrations may be improved if sub-lake algorithms are 
implemented. 

2. Little is known concerning the spatial and temporal distribu- 
tions of dissolved organic substance in Great Lakes waters. 
Since it is known that these substances can influence the 
upwelling spectra, measurement of dissolved organics should be 
included in any subsequent investigations which contain a 
remote sensing element. 

3. In the present program the sea truth sampling scheme consisted 
of detailed sampling at selected points. The scarcity of sea 
truth values limited our ability to place statistical confi- 
dence on much of the results. Future investigations of this 
type should look at transect and/or array collection schemes. 

4. While CZCS appears suitable to monitor whole lake conditions 
uncertainty remains as to how the CZCS will sample complex 
coastal waters. The density of supporting sea truth is related 
to the areal complexity of surface concentrations and the need 
to describe their spatial relationships. It is difficult, if 
not impossible, to directly sample complex waters and geograph- 
ically relate those features to coarse satellite observations. 
For purposes of understanding more completely how CZCS is 
averaging these features a study should be made where a high 
resolution remote sensor (e.g., Landsat-D TM) is used as an 
immediate stage for sampling. 
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INHERENT WATER OPTICAL PARAMETERS 
AS MEASURED BY THE NASA LANGLEY 
MOBILE WATER OPTICAL LABORATORY 


This portion of Appendix A contains the NASA Langley Research 
Center measured optical parameters for eighteen samples as Tables A1 
through A28 [10]. Measured optical parameters include: 

1. Volume scattering function (theta) 

2. Beam attenuation (alpha) 

3. Absorption (A) 

These measurements were made at eight wavelengths in 50 m increments 
from 450 nm through 800 nm. The above measurements were made with three 
separate instruments as described by Whitlock [31]. Three additional 
samples were collected in Lake Erie for optical analysis, but contained 
sufficient sediment to make measurement unreliable. Results from Tables 
A1 through A9, All, and A18 were used in the regressionn analysis 
results reported in section 6.1. 
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TABLE A1 



13 July 1980 Lake Erie Rocky River Mouth 


theta 

4S0. 

500. 

550. 

600. 

650. 

700. 

750. 

800. 

.374 

. 10631E+04 

.924366*03 

.801346*03 

. 74986E+03 

.694746*03 

.684676*03 

.63945E*03 

.600286*03 

. 7S1 

. 348646*03 

.301706*03 

.272316*03 

.25355E*03 

.236466*03 

.22895E*03 

.205086*03 

. 196096*03 

1.490 

. 944606*02 

. 94242E *02 

.353886*02 

.588756*02 

.786496*02 

•72760E*02 

.689026*02 

.630616*02 

as. eee 

. 659786*00 

.657796*00 

.626386*00 

.577246*00 

.533436*00 

.540246*00 

.495136*00 

.457416*00 

30.000 

.424286*00 

. 41434E*00 

.376766*00 

.369026*00 

.332916*00 

.344416*00 

.308946*00 

. 286686*00 

45 . 000 

. 14495E*00 

.140666*00 

.128046*00 

.140296*00 

.111296*00 

• 1 17S36*00 

.102626*00 

. 100166*00 

60 . 000 

.654606-01 

.625036-01 

•59830E-01 

.S4660E-01 

.517466-01 

. 52047E-01 

.483556-01 

. 43204E-01 

75.000 

.343186-01 

.326326-01 

.31487E-01 

.156666-01 

.268076-01 

. 28042E-01 

.255596-01 

.230056-01 

90.000 

.222466-01 

.218146-01 

. 19630E-01 

.212486-01 

.167376-01 

.172276-01 

. 15219E-01 

.144676-01 

105.000 

. 17039E-01 

. 167466-01 

. 15439E-01 

. 168506-01 

. 12921E-01 

. 13677E-01 

.116016-01 

.115696-01 

150.000 

.156736-01 

.149686-01 

. 13668E -01 

.148346-01 

. 1 1607E-0 1 

.120076-01 

.101516-01 

. 107106-01 

1 35 . 000 

.154306-01 

.150676-01 

. 13936E-01 

. 146S9E-01 

.117826-01 

. 128146-01 

. 104676-01 

.103796-01 

150.000 

.185416-01 

.177746-01 

. 16923E-01 

.149056-01 

.139576-01 

.144306-01 

•121076-01 

. 124656-01 

1S5.000 

.205676-01 

. 199396-01 

. 179806-01 

.170086-01 

. 15842E-01 

.167806-01 

.133966-01 

.133366-01 

ALPHA 

3.97 

3.39 

3.07 

2.93 

2.76 

2.83 

4.84 

4.14 

A 

.66 

.42 

.41 

.45 

.44 

.64 

2.60 

1 .95 


TABLE A2 

14 July 80 Lake Erie West Basin Hydra I, Sample 1 


TH6TA 

450. 

500. 

550. 

600. 

650. 

700. 

750. 

800. 

.374 

.842866*03 

.689166*03 

.575976*03 

.488406*03 

.43846E*03 

.409746*03 

.369636*03 

.330256*03 

•751 

.392916*03 

.329396*03 

.287936*03 

.241816*03 

.21605E*03 

. 199506*03 

. 17695E*03 

.160406*03 

1.490 

.141266*03 

.137776*03 

. 12052E*03 

.793636*02 

. 101 45E + 03 

.883716*02 

.794886*02 

.716056*02 

25.000 

.951936*00 

.95332E*00 

.864176*00 

.807146*00 

. 72726E + 00 

. 72427E +00 

•65819E+0O 

.586986*00 

30.000 

.591656*00 

.594756*00 

.530756*00 

. 5041 8E *00 

,44818E*00 

. 44 34 1 6 *00 

. 40249E *00 

.382746*00 

45.000 

.195286*00 

.199236*00 

. 17483E*00 

.167756*00 

. 15160E*00 

. 15359E+00 

. 13904E*00 

.130206*00 

60.000 

.869116-01 

.85762E-01 

.77916E-01 

.728616-01 

. 65010E-01 

.67189E-01 

.603236-01 

.542916-01 

75.000 

.436976-01 

.430906-01 

. 40874E-01 

.428116-01 

.342096-01 

.350026-01 

.313346-01 

.292646-01 

9e.000 

.282026-01 

.279126-01 

.25484E-01 

.267586-01 

.210476-01 

.216926-01 

.190386-01 

.190496-01 

105.000 

.216916-01 

.214356-01 

.185506-01 

.203516-01 

. 160996-01 

.163296-01 

. 1 45296-01 

.144016-01 

120.000 

.193746-01 

.190646-01 

.164276-01 

. 185236-01 

.143106-01 

.143806-01 

.129016-01 

. 128756-01 

135.000 

.195806-01 

.188816-01 

.163446-01 

. 18903E-01 

. 14046E-01 

.149146-01 

. 13340E-01 

. 12364E-01 

150.000 

.236916-01 

•22212E-01 

.20O32E-01 

. 18872E-01 

.165906-01 

. 177526-01 

. 155846-01 

. 1 476 4 E -0 1 

155.000 

.26010E-01 

.241176-01 

.21040E-01 

.275276-01 

.187696-01 

.204696-01 

.175476-01 

. 18565E-01 

ALPHA 

5.31 

4.67 

4.10 

3.78 

3.51 

3.51 

5.36 

4 . 49 

A 

.87 

.62 

.55 

.59 

.55 

.71 

2.65 

£.10 
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TABLE A3 

14 July 80 Lake Erie West Basin Hydra I, Sample 2 


theta «se. see. sse. 6ee. ese. ?ee. ?se. see. 

.374 .136015*04 . 1 1673E+64 .103625*04 .89742E*03 .815625*03 .80?12E*03 .730275*03 .678795*03 

.751 .508365*03 .437615*03 .38099E*03 .33016E*03 .29961E*03 .293825*03 .26456E*03 .241945*83 

1.498 .157285*83 .163235*83 .14568E+83 .982245*82 .123545*83 .113975*83 .181525*83 .915965*82 

25.888 .719425*88 .71428E*88 .65451E*00 .619815*88 .536395*88 .56257E*88 .49831E+00 .42984E*88 

38.888 .461885*88 .45154E*88 .482255*88 .38831E+88 .329645*88 .334985*88 .318875*88 .265485*88 

45.888 .168375*88 .167225*88 .143625*88 .158285*88 .114185*88 .122375*88 .116375*88 .183565*88 

68.880 .712S2E-81 .739535-01 .63920E-81 .59538E-81 .538835-01 .545045-01 .472665-01 .435085-01 

75.000 .376775-81 .481895-01 .333135-81 .40975E-01 .27836E-01 .382625-01 .267985-81 .254035-01 

90.000 .246035-01 .249165-01 .214635-01 .22268E-01 .18358E-01 .19436E-01 .172355-01 .156315-01 

105.000 .187455-01 .197215-01 .16G03E-01 .19989E-01 .14065E-01 .1S6S6E-01 .12435E-01 .121825-01 

120.800 .178895-01 .177245-01 .1S3415-01 .16616E-01 .122595-01 .133855-01 .120335-01 .183745-01 

135.000 .179415-81 .185065-01 .147825-01 .165505-01 .13292E-01 .132895-01 .138105-81 .110*55-01 

150.000 .225115-01 . 22627E-01 .186825-01 .183485-01 .16474E-01 .16807E-01 .15313E-01 .145165-01 

155.000 .25406E-01 .258095-01 .215495-01 .259985-01 -184535-01 .194175-01 .16700E-01 .165385-01 

ALPHA 4.95 4.29 3.69 3.55 3.37 3.37 5.27 4.30 

A .73 .51 .44 .53 .45 .68 2.66 2.08 


TABLE A4 

15 July 80 Lake Erie Cattaraugus Creek Mouth, Sample 1 


THETA 450. 580. 550. 608. 650. 700. 750. 800. 

.374 . 149845*04 .118165*84 .97840E*03 .8447^e*83 .757145*83 .768425*03 .706025*03 .646145*03 

.751 .691865*03 .576145*83 .500815*03 .441825*03 .412825*03 .4C498E*03 .367365*03 .347635*83 

1.498 .223765*03 .215195*03 .189635*83 .133775*03 .16609E+83 .16000E*03 -14643E*03 .136735 >83 

25.888 .191955*81 .191375*01 .190315*01 .173175*01 .174515*01 .187935*01 .166255*01 .16938E*01 

30.080 .122215*01 . 12402E+01 .12056E+01 .18577E+01 -11209E+01 .113115*01 .10443E*01 .102645*01 

45.000 .43221E+88 .455585*08 .436695*08 .421155*08 .39521E+00 .44337E*00 .361335*00 .369148*08 

60.000 .206405*00 .199505*00 .200335*00 .17805E*00 .179965*00 .185125*08 .159905*00 .153075*00 

75.000 .106345*88 .110695*00 .102845*00 .10266E*00 .91528E-01 .946585-01 .882Q1E-01 .823755-01 

90.000 .697245-81 . 70821E-81 .676395-01 .6S192E-01 .573785-01 .607935-01 .537485-01 .502915 01 

105.000 .531855-01 .529925-81 .514875-01 .470805-01 .435015 01 .451005-01 .484565-01 .406935-01 

120.000 .470925-01 .457795-01 .43054E-01 .42612E-01 .38020E-01 .986895-01 .940O7E-01 -34623E 01 

135.000 .456535-01 .434405-81 .428885-01 .408765-01 .377285-01 .367065-01 .32324E-01 .349065-01 

150.000 .522325-01 .50008E-01 .465165-01 .391625-01 .425015-01 .43251E-01 .37441E-01 .384135-81 

155.000 .594375-01 .560445-01 .51522E-01 .42S41E-01 .48563E-01 .453095-01 .41743E-01 .428945-01 

ALPHA 9.76 9.11 7.96 8.07 7.54 ^.19 9.10 8.29 

A 1.40 1.12 1.04 1.07 1.03 1.22 3.04 2.47 
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TABLE A5 

15 July 80 Lake Erie Cattaraugus Creek, Sample 2 


THETA 

450. 

500. 

550. 

800. 

650. 

700. 

750. 

800. 

.374 

. 10341E*04 

.8S215E*03 

.724096*03 

•62661E+03 

.S5149E+03 

. 53595E+03 

.504306*03 

. 50670E *03 

.751 

•39347E*03 

. 338726*03 

. 30097E*03 

.274616*03 

.253418*03 

. 24058E*03 

.222088*03 

. 2 1 347E *03 

1.490 

. 10907E*03 

. 10765E*03 

.963226*02 

. 64644E *02 

. 87489E*02 

.82990E«02 

. 77824E+02 

. 74965E *02 

2 s. 0 oo 

.878636*00 

. 903S6E+00 

. 857006*00 

. 81 420E *00 

. 75353E*00 

.77364E+00 

. 74 1 74E*00 

.67473E*00 

30.600 

.S8427E*00 

.S4S66E+00 

. 525566*00 

.511 59E*00 

.487956*00 

. 49214E+00 

. 46866E+00 

. 423996*00 

4S.O00 

. 205746*00 

.20420E400 

.188656*00 

.20696E+00 

. 16480E+00 

. 1741 7E*00 

. 15751E*00 

. 14392E*O0 

£0.000 

•96029E-01 

.95991E-01 

.87501E-01 

.84970E-01 

•77167E-01 

.777806-01 

. 72054E-01 

. 70046E-0 1 

75.000 

.51807E-O1 

•48915E-01 

.50300E-01 

•50312E-01 

.40725E-01 

.40564E-01 

. 38092E-01 

•3S441E-01 

90.000 

. 32572E-01 

. 33086E-01 

•27919E-01 

. 33676E-01 

.25476E-01 

. 26393E-01 

.24344E-01 

•22949E-0! 

105.000 

.25320E-01 

.26050E-01 

.22258E-01 

. 2381 16-01 

•20463E-01 

. 19842E-01 

•201S4E-01 

. 18428E 01 

120.000 

. 22231E-01 

.21G07E-01 

. 20585E-01 

. 22069E -01 

. 1G273E-01 

. 1784 7E -01 

. 16019E-01 

. 14242E-01 

135.000 

•21757E-61 

.20940E-01 

•21704E-01 

. 20619E-01 

. 16914E-01 

. 17266E-01 

. 16534E-01 

. 15729E-01 

150.000 

•24105E-01 

.2491 IE-01 

.22662E-01 

•2126OE-01 

. 1951 IE-01 

. 1996 r -E-01 

. 18408F-01 

. 19737E-01 

155.000 

.28081E-01 

.27240E-01 

. 23818E-01 

. 26432E-01 

•22177E-01 

. 23229E-01 

. 2 1059E -0 1 

•20439E -01 

ALPHA 

4.03 

3.38 

2.93 

2.81 

2.73 

2 . 96 

5.03 

4 .02 

A 

.79 

.70 

.71 

.68 

.64 

.85 

2.83 

2.17 


TABLE A6 

16 July 80 Lake Erie, North Shore, Sample 1 


THETA 

450. 

500. 

550. 

600. 

650. 

700. 

750. 

800. 

.374 

.629086*03 

.504866*03 

. 423S5E+03 

. 368876*03 

•32S45E*03 

.31823E*03 

.282406*03 .271896*03 

.751 

.334S4E+03 

.271846*03 

.229666*03 

.200416*03 

.180626*03 

.169216*03 

.153556*03 .141846*03 

1.490 

.13819E*03 

. 127S5E+03 

.108656*03 

.696346*02 

.887216*02 

.80420E*02 

.731446*02 .670006*02 

25.000 

.179686*01 

.174616*01 

. 16733E*01 

.138346*01 

.146476*01 

. 14795E*01 

.130386*01 . 11944E*01 

30.000 

.114586*01 

. 10773E*01 

.104346*01 

.851506*00 

.90749E*00 

.893286*00 

.831316*00 .7471SE*00 

45.000 

.38841E*00 

.361126*00 

.35207E+00 

.325296*00 

.304206*00 

.301476*00 

.281556*00 .2S8S9E*00 

60.000 

. 17347E+00 

. 16816E+00 

.146106*00 

.123156*00 

. 12854E+00 

.136526*00 

. 1 1898E*00 .111346*00 

7S.000 

•87231E-01 

.84438E-01 

.842696-01 

.809246-01 

•67471E-01 

.732536-01 

.610736-01 .584746-01 

90.000 

•59S33E-01 

•55337E-01 

.514496-01 

. 49589E-01 

. 4290SE-01 

.428176-01 

.389406-01 .370476-01 

105.000 

.43983E-01 

•39965E-01 

•37426E-01 

. 36807E-01 

. 33296E-01 

.34136E-01 

•29226E-01 .27673E-01 

120.000 

. 40227E-01 

•36576E-01 

.34103E-01 

. 34200E-01 

.287566-01 

•29754E-01 

.251116-01 .25020E-01 

135.000 

•37560E-01 

.367826-01 

.33000E-01 

.321626-01 

.282396-01 

.298856-01 

.23874E-01 .242886-01 

150.000 

. 43226E-01 

.398146-01 

.36415E-01 

.301236-01 

.306136-01 

.315406-01 

.271516-01 

.270586-01 

155.000 

. 46819E-01 

.425846-01 

.370826-01 

. 32538E-01 

. 32S94E-01 

.317106-01 

.27961E-01 .287706-01 

ALPHA 

7.40 

6.71 

5.60 

5.23 

4.89 

4.79 

6.59 

5.71 

A 

1.10 

.86 

00 

00 

.93 

.85 

1.03 

2.93 

2.22 
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TABLE A7 

16 July 80 Lake Erie, North Shore, Sample 2 



THETA 450. SM. 

.374 .79786E*03 .6407ge*03 

-751 . 489446*03 .396976+03 

1.4M .20327E+03 . 193666*03 

as.eee .a7498£+ei .esssec+oi 

30.000 . 176S3E*0t .168236+01 

45.000 .S9017E+00 .546996+00 

60.000 .264S4E+00 .249276+00 

75 . e«e . i3S94t*ee .i3iiiE+ee 

9e.eee .897i7E-ei .80S42E-ei 

les.eee .66744E-ei .6242SE-ei 

120.000 .S9i976-ei .ssieeE-ei 

i3S.eee .5?845E-ei . 522216-01 

tse.eee .664376-ei .S83i3£-ei 

155 . eee .?2669E-ei .654806-ei 

ALPHA 11.95 10.58 

A 1.41 1.22 


550. 600. 650. 

.528846+03 .457086+03 .408446*03 
.3363SE*03 .288636*03 .298706*03 
.164666*03 . 1088SE*03 .133336*03 
.245836*01 . 20395E*01 .20921E+01 
. 15291E+01 .127676*01 . 131036*01 
.507136*00 .462906*00 .463756*00 
.228626*00 .188926*00 .196616*00 
.111986*00 . 1 1230E +00 . 99788E-01 
.77900E-01 .70351E-01 .65744E-01 
•54330E-01 .52897E-01 .47669E-01 
.50981E-01 .49947E-01 .41U3E-01 
.48471E-01 . 4732SE-01 .40664E-01 
. 54039E-01 . 46616E-01 .46270E-01 
.56772E-01 . 48857E-01 .50060E-01 

9.12 8.27 7.62 

1.15 1.21 1.12 


700. 750. 800. 

•385036*03 .339936*03 .320086*03 
.243586*03 .215536*03 .194686*03 
.125766*03 .112386*03 .10122E*03 
.206416*01 . 183176+01 .171056*01 
. 12757E+01 .117256*01 .107266*01 
.438226*00 .38374E+00 .372576*00 
.199706*00 .171046*00 .165816*00 
.10593E*0e .890196-01 .843756-01 
. 64873E-01 .554256-01 .53191E-01 
• 48786E-01 .426986-01 .433666-01 
.441926-01 .37181E-01 .35933E-01 
. 42309E-01 .37617E-01 .360886-01 
.483726-01 .431246-01 .4160SE-01 
•S1954E-01 . 47734E-01 .44415E-01 

7.16 8.77 7.62 

1.36 3.28 2.56 


TABLE A8 

17 July 80 Lake Erie, West Basin, Hydra II, Sample 1 

theta 4S0. 500 . 550 . 600 . 6se. 700 . 7S0. 000 . 


.374 .10CS3E*04 .936706*03 

.751 .624306*03 .547406*03 

1.490 .275236*03 .278996*03 

25.000 .523136*01 .478526*01 

30.000 .335936*01 .29618E*91 

45.000 .118076*01 .111426*01 

80.000 .544996*00 .500126*00 

75.000 .285576*00 .264256*00 

90.000 .181346*00 .167326*00 

105.000 .137326*00 .127526*00 

120.000 .123296*00 .110846*00 

135.000 . 121516*00 .111156*00 

150.000 .149066*00 .134966*00 

155.000 . 175136*00 .154296*00 

ALPHA 25.68 23.21 

A 2.68 2.12 


.933986*03 . 786606*03 . 726126*03 
.491376*03 .4*5296*03 .415196*03 
.247906*03 .177576*03 .218846*03 
.510526*01 . 570086*01 .435526*01 
.310106*01 .366256*01 .277656*01 
.112606*01 .141266*01 .986046*00 
.497596*00 .555526*00 .4*0156*00 
.261966*00 .324966*00 .221886*00 
.169306*00 .204306*00 .135696*00 
.125736*00 .154466*00 .103646*00 
.112416*00 .138126*00 .953736-01 
.111146*00 .137206*00 .928556-01 
.134516*00 .134676*00 .116746*00 
. 1SS51E*00 .154266*00 .142796*00 

20.84 18.68 17.04 
1.79 1.83 1.75 


.7255tE*03 .657406*03 .643376*03 
.4] 13££«83 .371196*03 .354706*03 
.209676*03 . 1891 4E*03 .177576*03 
. 45593E*01 .39521E*01 .364136*01 
.289766*01 .251916*01 .243016*01 
.102466*01 .906186*00 .807686*00 
.453056*00 .383926*00 -36274E*00 
.242176*00 .202216*00 .190866*00 
.149746*00 .127*56*00 .11761E*00 
.112026*00 .94545E-01 .B8172E-01 
.991166-01 . 83933E-01 .780166-01 
.161656*00 .85*356-01 .823606-01 
. 121 1*E*00 . 10518E*00 .95*936-01 
. l*48SE*00 .123516*00 .11718E+00 

15.13 16.03 14.51 
1.94 3.45 2.87 
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TABLE A9 

17 July 80 Lake Erie, West Basin, Hydra II, Sample 2 


450. 500. 559 . 699. 656. 799 . 759 . 899. 

.92783E*03 .78S50E*03 .684196*03 .60S34E*03 .541066*03 .516036*03 .4S7S9€*03 .413836*03 

•41397E*03 .34273003 .29610E*03 .2S248E*03 .225966*03 .211286*03 . 191636*03 .174076*03 

. 14607E*03 . 141576*03 .123186*03 .802556*02 . 101306*03 .91S48E*02 .82S66E*02 .729986*02 

.888166*00 . 91494E *99 .833806*00 .77860E+00 .695736*00 .698826*00 .617516*00 .619036*00 

.576386*00 .577136*00 .520826*00 .48293E*00 .446656*00 .444S0E*00 .394206*00 .3S606E*00 

. 191626*00 . 194136*00 .17790000 .167246*00 .153896*00 . 1S644E+00 .136S2E*00 .129066*00 

.955506-01 .86S88E-01 .78608E-01 .738546-01 .68280E-01 .6GS86E-01 .589S4E-01 .53961E-01 

. 426306-01 .453866-01 .41743E-01 .44343E-01 .34867E-01 .353736-01 .32243E-01 -29S73E-01 

.27703E-01 .29336E-01 .2SS11E-01 .276496-01 .21861E-0I .222746-01 .19530E-01 .185596-01 

.226166-01 .21944E-01 .19443E-01 .20894E-01 . 16315E-01 . 171S3E-01 .14886E-01 . 1372SE-01 

. 1972GE-01 .20018E-01 .17S76E-01 .191S4E-01 .147S8E-01 .1SG30E-01 .13184E-01 .12701E-01 

. 19705E-01 .20214E-01 .17891E-01 .19970E-01 .14769E-01 .1S299E-01 .13789E-01 .12616E-01 

. 2401 IE-01 .23G94E-01 .2172SE-01 .20301E-01 .17324E-01 .190S7E-01 .162S9E-01 .1S780E-01 

.28724E-01 .27733E-01 .2S182E-01 .29308E-01 .21067E-01 .228G6E-01 .20S16E-01 .17193E-01 

4.93 4.28 3.76 3.4S 3.34 3.28 S.13 4.33 

.77 .SI .52 .52 .48 .64 2.70 1.97 


TABLE A10 

21 July 80 Lake Michigan, Green Bay 


450 . 500 . 550 . 600 . 650 . 700 . 750 . 800. 

.11468004 .93445043 .978S06*03 .87821E*03 .8S903E*03 .791986*03 . 832656*03 .733166*03 

.413136*03 .369S4E*03 .33265E*03 .29079E+03 .27120€*03 .25932E*03 .247976*03 .241096*03 

.100056*03 .103686*03 .911706*02 .G0298E*02 .850436*02 .77635E*02 .746306*02 .43707E+02 

.399766*00 .39854E*O0 .364546*00 .36827E*00 .306566*00 .329616*00 .29311E*00 .288346*00 

.247606*00 .252766*00 .234306*00 .215496*00 .201076*00 .213726*00 .18255E*00 . 17990E+00 

.838216-01 . 87207E-01 .821936-01 .862406-01 .701806-01 .72275E-01 .739376-01 .68306E-01 

. 37549E-01 . 39S78E-01 .355356-01 .31840E-01 .30928E-01 .382196-01 .31412E-01 .325116-01 

. 19191E-01 .20978E-01 .199436-01 .220146-01 .18405E-01 .20S44E-01 .18245E-01 .16320E-01 

.137776-01 .146126-01 .13403E-01 .138686-01 .108746-01 .122636-01 .108086-01 .10172E-01 

.114806-01 .101726-01 . 1 1065E-01 .988836-02 .784356-02 .822856-02 .853026-02 .101016-01 

.96126E-02 .933886-02 .948G8E-02 .109436-01 .717586-02 .803136-02 .654826-02 .720746-02 

.939166-02 .986286-02 .868986-02 .888906-02 .774846-02 .988666-02 .67051E-02 .796416-02 

. 1 1764E-01 . 1 1786E-01 .108436-01 .10655E-01 .935076-02 .10094E-01 .978996-02 .10S55E-01 

.139316-01 .145026-01 .137176-01 .11955E-01 .134266-01 .130386-01 .11645E-01 .11622E-01 

3.93 3.42 2.97 2.80 2.68 2.71 4.84 3.88 

.71 .33 .29 .29 .26 .41 2.38 1.81 


THETA 

.374 

.751 

1.490 

25.000 

30.000 

45.000 

60.000 

75.000 

90.000 

105.000 

120.000 

135.000 

150.000 

155.000 

ALPHA 

A 
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TABLE All 

21 July 80 Lake Michigan, Grand Haven, Sample 1 


THETA 

450. 

500. 

550. 

600. 

650. 

700. 

750. 

800. 

.374 

.751 

1.490 

25.000 

39.000 

45.000 

60.000 

75.000 

90.000 

105.000 

120.000 

135.000 

150.000 

155.000 

.28487E404 . 2354 8E 404 
.105676*04 .961386*03 
.273536403 .276566403 
.160466401 .177986401 
.102076401 .108596401 
.331406*00 . 34866E+00 
.142246*00 .155106*00 
•75335E-01 .823946-01 
.484976-01 .51567E-01 
.361846-01 .388136-01 
.328926-01 .348186-01 
.322396-01 . 34570E-01 
. 42173E-01 . 46364E-01 
.517896-01 .548266-01 

.205966*04 .190506*04 
.882386*03 .856806*03 
.253616*93 .18787E*03 
.177676*01 .167956*01 
.109546*01 .104S1E*01 
.3S147E*00 . 38413E*09 
. 15557E*00 . 15053E*00 
.795966-01 .874S3E-01 
.495246-01 . 56076E-01 
. 37059E-01 .410426-01 
. 32829E-01 . 3691 IE-01 
. 33302E-01 . 35892E-01 
. 430S9E-01 . 41749E-01 
. 53593E-01 .6559SE-01 

.184736*04 
.840026*03 
.241 446*03 
. 15500E*01 
.93256E*00 
. 309326*00 
. 134106*00 
.694116-01 
.420716-01 
. 31978E-01 
, 27769E-01 
.286446-01 
.367306-01 
.440436-01 

.172126*04 
.793936*03 
.24046E*03 
.165266*01 
.101056*01 
.344276*00 
.1471 1E*00 
.778856-01 
.461576-01 
.344086-01 
.304276-01 
.315726-01 
.421816-01 
. 52846E-01 

.156956*04 
.759776*03 
.2 39206*03 
.143196*01 
.9148S6*00 
.312646*00 
.130336*00 
.685786-01 
.417196-01 
.309236-01 
.270336-01 
.275996-01 
.376766-01 
.485546-01 

.144076*04 
. 7S319E*03 
■ 24330E *03 
.130166*01 
.824336*00 
.28007E+00 
.119476*00 
.611846-01 
. 37285E-01 
.282686-01 
.246476-01 
.24545E-01 
. 36356E -01 
. 46774E-01 

ALPHA 

13.66 

12.33 

10.78 

10.00 

9.35 

9.08 

10.90 

9.91 

A 

2.89 

2.06 

1.46 

1.37 

1 .35 

1.47 

2.94 

2.37 


TABLE A12 

22 July 80 Lake Michigan, Grand Haven, Sample 2 


THETA 45®. S0®. 550. 60®. 65®. 700. 75®. 80®. 


.374 . 761S5£*93 .639126*03 .55004E+03 .50416E*03 .459506*03 .445266*03 .429656*03 .40854E+03 

.751 .279006*03 .242696*03 .218086*03 .199876*03 .183316*03 .175756*03 .16S37E*03 .154136*03 

1.49® .735136*02 .74S61E*02 .67882E*02 . 45988E*02 .652136*02 .581506*02 .562776*02 .520246*02 

25.090 .267876*00 .264406*00 .24653E*00 .239816*00 .217046*00 .228956*00 .20132E*00 . 197[7F*00 

30.090 .170S2E*00 . 16185E*00 .1S344E+00 .14070E+00 .12782E*®0 .13199E+00 .118036*00 .11282E*00 

45.000 .51244E-01 .S2323E-01 .46725E-01 .54094E-01 .427S1E-01 .44631E-01 .376636-01 • 381- ir H--M 

60.000 .21574E-01 .21659E-01 .20003E-0J . 19511E-01 .17735E-01 .19083E-01 .163896-01 .14973!: 91 

75.000 . 112496-01 .111106-01 .10492E-01 .11730E-01 .902986-02 .944P5E-02 .79551E-C2 -79363E-02 

90.000 .77S95E-02 .730S6E-02 .65056E-02 .69956E-02 .56550E-02 .S8358E-02 .50131E-O2 .47326E-02 

105.000 .59337E -02 .56040E-02 .49427E-02 .52413E-02 .41974E-02 .42012E-02 .38490E-02 . 328446-02 

120.000 .54526E-02 .50345E-02 .44021E-02 .49887E-02 .429786-02 .38935E-02 .35849E-02 .34426E02 

135.000 .55300E-02 .52468E-02 .48388E-02 .51740E-02 .40648E-02 .41894E-02 .36529E-02 .36098E-02 

150.000 .815576-02 .77926E-02 .73216E-02 .68568E-02 .G6714E-02 .64169E-02 .58475E-02 .5S35GE-02 

155.000 . 10660E-01 . 10317E-01 .98170E-02 .11345E-01 .8256SE-02 .8S604E-02 .738386-02 .777966-02 

ALPHA 2.37 2.01 1.83 1.80 1.90 1.94 3.98 3.58 

A .11 .07 .11 .22 .15 .38 2.38 1.74 
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TABLE A13 

22 July 80 Lake Michigan, Grand Haven, Sample 3 


THETA 4SO. 5M. 

.374 . 471 44E+03 .396086+03 

.751 . 174516+03 . 152996+03 

1.490 . 42B69E+02 .443596+02 

25.000 . 141206+00 . 135676+00 

30.000 .84S76E-01 .804176-01 

45.000 .26736E-01 .27863E-01 

60.000 . 12041E-01 . 10980E-01 

75.000 .631866-02 .60905E-02 

90.000 .402036-02 . 38361E-02 

105.000 .33538E-02 .29987E-02 

120.000 .354706-02 .29192E-02 

135.000 . 37215E-02 .32907E-02 

150.000 .57833E-02 .S1445E-02 

155.000 .73600E-02 .68552E-02 

ALPHA 1.29 1.10 

A .01 .01 


550. 600. 650. 

.344926+03 .317326+03 .279666 *03 
.138886+03 .128306+03 .121936+03 
.403766+02 .26221E+02 .422706+02 
.123216+00 .107356+00 .108676+00 
.800386-01 .648626-01 .67552E-01 
.244026-01 .21678E-01 .22586E-01 
. 1081 IE-01 .910096-02 .91188E-02 
. 57238E-02 .49583E-02 .467746-02 
.363396-02 .33819E-02 -299126-02 
.26960E-02 .24261E-02 .25402E-02 
.26971E-02 .238546-02 .233636-02 
. 30 191 E -02 . 26705E-02 -23981E-02 
.479856-02 .36875E-02 .38261E-02 
.65961E-02 .63296E-02 .54646E-02 

.98 1.07 1.11 

.03 .11 .07 


700. 750. 800. 

.271036+03 .249216+03 .25640E+03 
. 1 1542E+03 .107736+03 .10175E+03 
.359026+02 .349756+02 .32791E+02 
.117496+00 .971536-01 .104786+00 
. 69284E-01 .602616-01 .60425E-01 
.227206-01 .19844E-01 .19025E-01 
.110896-01 .787466-02 .98865E-02 
.508656-02 .427926-02 .476966-02 
.292466-02 .259936-02 .28799E-02 
.227286-02 .199976-02 .21144E-02 
.218246-02 . 17628E-02 .210416-02 
. 25992E -02 .232826-02 .23929E-02 
.383896-02 .34068E-02 .386746-02 
•55489E-02 .50947E-02 .539686-02 

1.31 3.36 2.37 

.31 2.33 1.77 


TABLE A14 

July 80 Lake Michigan, Grand Haven, Sample 4 


THETA 450. S00. 550. 600. 650. 700. 750- 800. 


.374 
• 751 

1.490 

25.000 

30.000 

45.000 

60.000 

75.000 

90 . 000 

105.000 

120.000 

135.000 

150.000 

155.000 


.590696+03 .460336+03 
. 19853E+03 .167606+03 
. 47608E+02 . 48884E +02 
.24039E+00 .24463E+00 
.146316+00 .147846+00 
. 4S103E-01 .462116-01 
. 18959E-01 . 18616E-01 
. 97918E-82 .966646-02 
•60440E-02 .60921E-02 
. 48575E-02 . 47550E-02 
.449666-02 . 44296E-02 
•50694E-02 .49688E-02 
. 70822E-02 .66554E-02 
.91054E-02 .888416-02 


.418886+03 .39652E+03 
. 16094E+03 . 15008E+03 
-44825E+02 .29303E+02 
.23214E+00 .181556+00 
.135376+00 .114136+00 
. 43762E-01 .408766-01 
. 17416E-01 .145096-01 
•91571E-02 .84257E-02 
.553146-02 .509826-02 
.412316-02 .38731E-02 
.384976-02 .370656-02 
.417366-02 .375736-02 
. 55838E-02 .364476-02 
.742726-02 .887736-02 


.352226+03 .325216+03 
.145936+03 .137676+03 
.445906+02 .394976+02 
.195056+00 .207856+00 
.119146+00 .123896+00 
.379436-01 .382+36-01 
.152996-01 .173826-01 
.760616-02 .824116-02 
.472746-02 .499+16-02 
.365696-02 .40932E-02 
. 33408E-02 . 35642E-02 
•35646E-02 .382216-02 
.523446-02 .54830E-02 
.723186-02 .757436-02 


.292486+03 .288006+03 
.122366+03 .115236+03 
. 40692E+02 .373796+02 
.180986+00 .169956+00 
.114536+00 .105436+00 
.362926-01 .345726-01 
.144586-01 .154976-01 
.73759E-02 .83961E-P2 
.464206-02 .559726-02 
.372236-02 .483556-02 
.341406-02 . 46756E-02 
.332136-02 .420626-02 
.472076-02 .511446-02 
.676716-02 .706316-02 


ALPHA 2.06 1.82 1-60 1.61 1.64 1.79 3.78 3.21 


A .03 .00 .03 -15 .07 .30 2.26 1.87 
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TABLE A15 

23 July 80 Lake Michigan, Grand Haven, Sample 5 


THETA 

459. 

599. 

550. 

609. 

650. 

700. 

750. 

800. 

.374 

.229276+04 .187986+94 

.168686+04 

.145226+04 

.134036+04 

. 128876+94 

.117796+04 

. 114676+04 

. 7S1 

.763326+03 . 681336+93 

.648236*03 

.627956+03 

.603316+03 

.590436+03 

.555466+03 

.546806+03 

1.499 

.182346+93 . 167316+93 

.172936+03 

.126136+03 

. 165446+03 

.165976+03 

. 16005E + 03 

. 157856+03 

25.099 

.999696+99 .98547E+99 

.977076+00 

.930706+00 

.838206+00 

.880156+90 

. 85099E+00 

.786986+00 

39. e«e 

.571376+99 .604746+99 

. 58990E+00 

.555796+00 

.512766+00 

.556336+00 

.516506+00 

.517046+00 

45.999 

. 178526+99 .197446+99 

. 191676 + 00 

.204286+00 

. 165696+00 

.175516+00 

. 16564E + 00 

. 15363E+00 

69.999 

•7793S6-91 .839466-91 

.795806-01 

. 78531E-01 

.699346-01 

. 74223E-01 

.69869E-01 

•69346E-01 

7S.999 

.495986-91 .436876-91 

.492896-01 

.448126-01 

.349416-01 

.386986-01 

.179406-01 

.328706-01 

99.999 

.247956-91 .264736-91 

.248776-01 

.265836-01 

.216266-01 

.235686-01 

.201436-01 

. 2O266E-01 

195.999 

.199196-91 .197896-91 

. 18240E-01 

.204196-01 

.159426-01 

. 165226-01 

. 15408E-01 

.156316-01 

129.999 

. 16672E-91 .176966-91 

.157406-01 

.182266-01 

.139116-01 

.146746-01 

.131306-01 

.134386-01 

135.999 

.177626-01 .184266-01 

.179706-01 

.190286-01 

.147256-01 

. 156816-01 

.142216-01 

.141046-01 

159.999 

.230336-91 .248396-01 

.229796-01 

. 22623E-01 

.196226-01 

. 20706E-0 1 

.180726-01 

. 167446-01 

155.999 

.322596-01 .313246-91 

.297046-01 

.381516-01 

.26430E-01 

.291026-01 

.253366-01 

.24933E-01 

ALPHA 

7.35 

6.75 

6.35 

5.95 

5.80 

5.74 

7.52 

6.79 

A 

1.46 

.97 

.75 

.76 

.73 

.88 

2.66 

2.16 


TABLE A16 

24 July 80 Lake Superior, Duluth, Sample 1 


THETA 

459. 

500. 

550. 

600. 

650. 

700. 

750. 

■_ 3 > > 

.374 

.709266+03 

.567426+03 

.462476+03 

.406716+03 

.360416+03 

.346236+03 

.30575E+03 

. rr ♦ o } 

.751 

.321876+03 

•26361E+03 

.222176+03 

. 19209E + 03 

. 174196+93 

. 16.7636 + 03 

. 14',?;jE + 03 

. 1 ■•_*:' -f + 03 

1.499 

. 1 1302E + 03 

.103626+03 

•87252E+02 

.570006+02 

. .’4 )616 + 08 

.60 P63E+02 

.612356+08 

.'66+02 

25.990 

.166266+01 

. 16342E+01 

. 144406 + 01 

-11647E+01 

. 1282S6+01 

. 1 15406+01 

. 9‘-0 1 06 + 00 

. 3 .r- r -..r + 00 

39.099 

. 1 1 369E + i4i 

. 104196+01 

. 99048E+00 

.790376+00 

. 81 485E+00 

.7705 IT +09 

.6433(16+00 

.6131-L+00 

45.099 

.432686+09 

.398416+00 

•3?67e6+O0 

.357386+00 

.314656+00 

.30236.6*00 

.258006+00 

. 24+10 1 f * 00 

60.900 

.197086+ . ' 

1 9 7 1 76 + 00 

. 18351E+00 

. 154566*00 

. 149306+00 

. 141106+00 

. 122951+00 

. 1 8 1 . 176 + 00 

75 . 90 u 

. 10705E+00 

. 1049SE+00 

. 100856+00 

.970156-01 

. 86201E-01 

.799556 01 

.711666-01 

.603636- 01 

99.009 

.720266-01 

.685596-01 

.65361E-01 

.652916-01 

. 55604E-01 

.539246-01 

. 461H9L 01 

- >88 : : 1 C 01 

195.000 

.551936-01 

.5523SE-01 

.532606-01 

. 501 29E -01 

. 424386-01 

.426536-01 

. 391006-01 

. }, .>‘1006 +4 1 

120.000 

. 50407E-01 

.490956-01 

.457336-01 

.446326-01 

. 39865E -01 

.383236-01 

. 3 1 4756-01 

.297776 01 

135.009 

.485846-01 

.483766-01 

.457646-01 

.434456-01 

•38554E-01 

•36560E-01 

. 307726-01 

.280196-01 

150.909 

.579286-91 

.556756-01 

.503306-01 

.422136-01 

. 42 1 91E-01 

. 40837E-O1 

. 349476-01 

. 9871 36 -01 

155.000 

.683346-01 

.608706-01 

.577536-01 

. 492426-01 

. 43352E-01 

. 461656-01 

. 39+3996 -01 

. "1308116 -0 1 

ALPHA 

7.98 

6.34 

5.14 

4.47 

4.09 

4.03 

5 . 79 

c ,.0 

A 

2.44 

1 .66 

1 . 31 

1.11 

1 .01 

1 . 05 

8.87 

2 . 4 J 
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TABLE A17 

24 July 80 Lake Superior, Duluth, Sample 2 


THETA 

.374 

.751 

1.499 

25. 999 

39.999 

45.999 

69.999 

75.999 

99.999 

195.999 

129.999 

135.999 

159.999 

155.999 

ALPHA 


459. 

.498256*03 
. 239596*93 
.794356*02 
. 101536*01 
.661566*00 
.2479SE*00 
.11299099 
.61S46E-91 
.409376-01 
.31644E-91 
.29928E-91 
•28732E-91 
. 34432E-91 
.41579E-91 

4.95 

1.23 


599. 

.49S26E*93 
. 1871 16*93 
•72S33E*02 
.986866*09 
.653936 +00 
.241686499 
. 1 16846*00 
.634686-91 
.425876-91 
.339916-91 
.295696-91 
.29791E-91 
.344126-91 
.498526-91 

4.14 

.81 


559. 

.313566493 

.1S367E493 

.S9921E492 

.884176499 

.58921E+99 

.22297E+99 

.19328E499 

.586696-91 

.385266-91 

.397446-91 

.273316-91 

.26794E-91 

.39917E-91 

.357546-91 

3.59 

.68 


699. 

.264396493 
.138926403 
. 36478E+92 
.78455E409 
.538476400 
.240546400 
.994496-01 
.638946-91 
.412636-01 
. 33098E-01 
.29772E-01 
■28714E-01 
.289226-01 
. 35510E-01 

3.18 

.65 


659. 

.22713E403 
.121886403 
.595906402 
.700846400 
.485416400 
. 184366400 
. 88046E-01 
.472346-01 
.311966-01 
.249906-01 
. 2280S6-01 
.220556-01 
.248336-91 
.308526-01 

2.95 

.55 


700. 

.212346403 
. 114746403 
.44637E402 
. 67888E+00 
.456356400 
.17565E400 
.830536-01 
.491616-01 
.398636-01 
.23865E-01 
.221586-01 
.21965E-01 
,253606-01 
.296416-01 

2.80 

.74 


750. 

.189096403 
.106666403 
.427696402 
.591956400 
.397006400 
. 154556*00 
.714196-01 
.412726-01 
.267476-01 
.20489E-01 
.185786-01 
. 184226-01 
.221766-01 
.265796-01 

4.68 

2.67 


890. 

.179326*03 

.952826*02 

.389006*92 

.531756*00 

.356426*00 

.137526*00 

.681656-01 

.355456-01 

.240656-01 

.194166-01 

.171276-01 

.174956-01 

.214846-01 

•2S120E-01 

3.91 

1.99 


TALBE A18 

25 July 80 Lake Michigan, Grand Haven, Sample 6 


THETA 459. 599. 550 . 699. 


659. 


790. 759. 809. 


.374 

.751 

1.499 

2S.999 


39.999 

45.999 

69.999 


75.999 


90.990 

195.999 

129.990 

135.990 
150.090 

155.990 


.267996*94 
.811706*03 
.199126*93 
.852226*99 
.527736*99 
.161516*99 
.697446-01 
.363186-91 
. 228696-01 
.179936-91 
.155486-01 
. 164696-91 
.21375E-01 
.285626-01 


.224616*94 
.758216*93 
.196466*93 
.967656*99 
.587236*99 
. 188926*99 
.894836-91 
.409196-01 
.246126-01 
. 1834SE-01 
. 16299E-01 
.168806-01 
.220696-01 
.270266-91 


.1928SE*94 
.73S31E+03 
.184116*03 
.94833E400 
.S7403E400 
.183616*09 
•75782E-01 
.3S288E-01 
•24123E-01 
.172516-01 
. 17576E-01 
.160166-01 
.213926-01 
.270696-91 


.163636*94 
.70413E403 
.134486*03 
.885156*00 
.53036E*09 
.193186*00 
.712236-01 
. 41853E-01 
.247556-01 
.182276-01 
.171026-01 
.172386-01 
.208266-01 
.35789E-01 


. 14481E*04 
.674646*03 
.175296*03 
.826316*00 
.506806*00 
.160396*90 
.672636-01 
.337916-01 
.201086-01 
.148296-01 
.131146-01 
. 13805E-01 
.177006-01 
.23180E-01 


.135146*04 
.64671E*93 
•17269E*03 
.909116*00 
.552306*00 
.17843E400 
.717956-01 
. 37217E-01 
.221666-01 
. 16668E-01 
.144896-01 
. 15177E-01 
.199976-01 
.25912E-01 


. 121376*04 
.628826*03 
.172436403 
.897906*90 
.50784E+00 
. 1608SE40O 
.67904E-01 
.334166-01 
.20398E-01 
. 14621E-01 
. 12680E-01 
.138056-01 
. 17439E-01 
.214096-01 


.121096*04 
.643786*93 
. 17736E*03 
.790166*00 
.50286E*00 
. 15557E+00 
.677686-01 
.33137E-01 
.196886-01 
.149436-01 
.127716-01 
.133946-01 
.173806-01 
.238176-01 


ALPHA 8.28 


7.50 6.90 6.50 


6.33 6.36 


8.27 7.58 


A 


1.80 1.17 .80 


.80 


.80 


.92 2.59 2.21 


129 


NASA SEA TRUTH SAMPLING DATA 


This portion of Appendix A contains sea truth sampling data as 
collected and analyzed by the University of Michigan Research Center for 
the Great Lakes [20]. Parameters measured include: 

Physical: Secchi depth, sample temperature, wave height, wind 
speed, wind direction, water color, sky color, and cloud cover. 

Biochemical: Chlorophyll pigments (A, A2, B, and C), Phaeophytins, 
filter residue (total, ashed, volatile), and particulate organics as 
carbon, hydrogen and nitrogen fractions. 


These data are reported in Tables A19 through A32. 
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APPENDIX B 

CZCS COLOR CODED RADIANCE, CHLOROPHYLL-A PIGMENT, 
SUSPENDED SEDIMENT, AND SURFACE TEMPERATURE MAPS 
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1. Spatial Comparison of Atmospheric Correction Algorithms for 
Lake Superior 

2. Spatial Comparison of Atmospheric Correction Algorithms for 
Lake Michigan 

3. Chlorophyll -a and Suspended Sediment Concentration Maps, Lake 
Huron Subscenes 

4. Chlorophyll-a and Suspended Sediment Concentration Maps, Lake 
Michigan Subscenes 

5. Chlorophyll-a and Suspended Sediment Concentration Maps, Lakes 
Erie and Ontario 

6. Chlorophyll-a and Suspended Sediment Concentration Maps, Lake 
Superior Subscenes 

7. Upwelling Radiance Map, Lake Huron at 443 and 550 nm 

8. CZCS Surface Temperature Maps, Lakes Michigan and Erie 

9. CZCS Surface Temperature Maps, Lakes Superior and Huron 
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Terim 


PLATE I 


SPATIAL COMPARISON OF ATMOSPHERIC CORRECTION 
ALGORITHMS FOR THE LAKE SUPERIOR SUBSCENE 

This plate shows a series of color coded radiance maps 
o 

(mw/cm .sr.p) for each of the first four CZCS channels as collected on 
June 11, 1980. Shown are the original uncorrected radiance image which 
does not include any adjustment for sensor deterioration, and the 
corrected radiance map for each of the following three algorithms: 
Smith and Wilson [5] (Algorithm 2), Jain's pseudo optical depth 
algorithm [6] (Algorithm 3), and the residual component algorithm 
(Algorithm 4). The corrected radiances are shown as transformed to a 
nadir view geometry and a Rayleigh atmosphere. 
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PLATE II 

SPATIAL COMPARISON OF ATMOSPHERIC CORRECTION 
ALGORITHMS FOR THE LAKE MICHIGAN SUBSCENE 

This plate shows a series of color coded radiance maps for each of 
the first four CZCS channels as collected on June 11, 1980. Shown are 
the original uncorrected radiance image which does not include any 
adjustment for sensor deterioration, and the corrected radiance map for 
each of the following three algorithms: Smith and Wilson [5] (Algorithm 

2), Jain's pseudo optical depth algorithm [6] (Algorithm 3), and the 
residual component algorithm (Algorithm 4). The corrected radiances are 
shown as transformed to a nadir view geometry and a Rayleigh Atmosphere. 
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PLATE III 


CHLOROPHYLL-A AND SUSPENDED SEDIMENT 
CONCENTRATION MAPS, LAKE HURON SUBSCENES 

o 

Color coded chlorophyll -a pigment concentration (C mg/m ) maps are 
shown in the left panel for Lake Huron (upper), Saginaw Bay (lower 
left), and Lake St. Clair (lower right). Corresponding suspended 
sediment maps (SM mg/ji) are shown in the right panel. These maps were 
created from the June 11, 1980 CZCS scene. The suspended sediment 
algorithm was parameterized to map the high concentrations of suspended 
sediments in Saginaw Bay. For the Saginaw Bay subscene the predicted 
very low concentration of less than 0.1 mg/i, SM outside the Bay is not 
considered representative of these waters and was a consequence of the 
algorithm parameterization. 
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PLATE IV 


CHLOROPHYLL-A AND SUSPENDED SEDIMENT 
CONCENTRATION MAPS, LAKE MICHIGAN SUBSCENES 

3 

Color coded chlorophyll -a pigment concentration (C mg/m ) maps are 
shown in the left panel for Green Bay (left), lower Lake Michigan 
(middle) and Lake Michigan (right). Corresponding suspended sediment 
maps (SM mg/t) are shown in the right panel. These maps were created 
from the June 11, 1980 CZCS scene. The suspended sediment algorithm was 
parameterized to map the high concentrations of suspended sediments in 
Green Bay. For the Green Bay subscene the predicted very low concentra- 
tion of <0.1 mg/£ SM outside the Bay is not considered representative of 
these waters and was a consequence of the algorithm parameterization. 
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PLATE V 


CHLOROPHYLL-A AND SUSPENDED SEDIMENT 
CONCENTRATION MAPS, LAKES ERIE AND ONTARIO 

•3 

Color coded chi orophyl 1 -a pigment concentration (C mg/m ) maps are 
shown in the left panel for Lake Erie (upper left), weestern Lake Erie 
(upper right), and Lake Ontario (lower middle). Corresponding suspended 
sediment maps (SM mg/s,) are shown in the right panel. These, maps were 
created from the July 14, 1980 CZCS scene. The suspended sediment 
algorithm was parameterized to map the high concentrations of suspended 
sediments in western Lake Erie. For the Lake Erie and western Lake Erie 
subscene the predicted very low concentration <0.1 mg/£ SM outside the 
western basin is not considered representative of these waters but were 
a consequence of the algorithm parameterization. In addition the 
central portion of Lake Ontario is cloud covered and the mapped 
concentrations are, therefore, not reliable for this area. 
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PLATE VI 


CHLOROPHYLL-A AND SUSPENDED SEDIMENT 
CONCENTRATION MAPS, LAKE SUPERIOR SUBSCENES 

Color coded chlorophyl 1-a pigment concentration (C mg/m ) maps are 
shown in the left panel for Lake Superior (upper). Black Bay (lower 
right), and western Lake Superior (lower left). Corresponding suspended 
sediment maps (SM mg/fc) are shown in the right panel. These maps were 
created from the August 6, 1980 CZCS scene. The suspended sediment 
algorithm was parameteri zed to map the high concentrations of suspended 
sediments in Black Bay. The predicted very low concentration of <0.1 
mg/£ SM outside the Bay is not considered representative of these waters 
and was a consequence of the algorithm parameterization. 
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PLATE VII 


UPWELLING WATER RADIANCE MAPS FOR 
LAKE HURON SUBSCENE AT 443 AND 550 NM 

O 

This plate shows a color coded water radiance map (mw/cm t sr . u) 
for CZCS bands 1 and 3 (443 and 550 nm). The 443 nm radiance maps are 
shown on the left panel and the 550 nm maps in the right panel. The 
upwelling surface radiances as calculated by the atmospheric correction 
processing have been normalized to a nadir view angle at the top of a 
Rayleigh atmosphere. The pixel to pixel radiances depict variations in 
the surface upwelling radiance. 
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PLATE VIII 


CZCS SURFACE TEMPERATURE MAPS FOR 
LAKES MICHIGAN AND ERIE SUBSCENES 

This plate shows surface temperature over the range from 11 to 25 
°C as 1 and 2 degree level slices. CZCS thermal band radiances were 
calibrated with the available surface measurements. A linear regression 
model was developed for these data and used to extrapolate to all 
observed radiance levels. Thus there is greater uncertainty in 
predicted temperatures at both ends of the scale. The thermal maps on 
the left are for Lake Erie (upper left), western Lake Erie (upper right) 
and Lake Ontario (lower middle). Note that the upper portion of Lake 
Ontario was partially cloud covered for this July 14, 1980 scene. The 
maps on the right show Green Bay (right), lower Lake Michigan (middle) 
and Lake Michigan (right). The thermal patterns in Lake Michigan have 
been shown for this June 11, 1980 scene to be highly correlated to lake 
bathymetry. 
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PLATE IX 


CZCS SURFCE TEMPERATURE MAPS FOR 
LAKES SUPERIOR AND HURON SUBSCENES 

This plate shows surface temperature in over the range from 7 to 21 
°C as 1 degree increments. CZCS thermal band radiances were calibrated 
with the available surface measurements. A linear regression model was 
developed for these data and used to extrapolate to all observed 
radiance levels. Thus there is greater uncertainty in predicted 
temperatures at both ends of the scale. The thermal maps on the left 
are for Lake Superior (upper), western Lake Superior (lower left) and 
Black Bay (lower right) and were made from the August 6, 1980 scene. 
The maps on the right show Lake Huron (upper), Saginaw Bay (lower left) 
and Lake St. Clair (lower right). Note that the upper right hand 
section of Lake Huron is partially cloud covered for this June 11, 1980 
scene. 
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